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ABSTRACT 

Composite  strength  is  an  inverse  function  of  the  size  of  the  com- 
posite. As  the  use  of  composites  expands  into  larger  applications, 
such  as  airplane  wings,  missile  components,  and  ship  superstructures, 
the  ability  to  accurately  predict  composite  performance  for  large 
applications  has  become  more  important.  The  composite  failure  pro- 
cess is  sequential  and  initiates  with  early  breaking  of  the  weak  fibers. 
Concentration  of  breakage  sites  accumulates  and  leads  to  ultimately 
catastrophic  failure.  Prestressing  fibers  prior  to  solidification  of  the 
matrix  has  been  demonstrated  to  increase  the  reliability  of  the  com- 
posite by  minimizing  the  spatial  concentration  of  the  breakage  sights. 

This  study  concentrates  on  quantifying  the  level  of  preload  and 
gauge  length  to  optimize  the  prestress  effect.  Computer  simulations  of 
graphite  bundle  tests  were  used  to  validate  the  data  analysis  method- 
ologies applied  to  actual  AS-4  graphite  bundle  tests.  The  actual 
experimental  results  are  consistent  with  computer-simulated 
behaviors. 
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I.     INTRODUCTION 

In  the  1980s,  the  capability  of  military  equipment  has  increased 
exponentially.  Modern  aircraft  fly  further  and  carry  greater  payloads, 
missiles  have  larger  warheads  and  greater  ranges,  tanks  have  bigger 
guns  and  tougher  armor,  and  ships  have  more  capable  but  lighter 
weapon  systems  along  with  improved  armor  plating.  Engineering 
improvements  in  the  areas  of  electronics  and  composite  materials 
have  been  on  the  leading  edge  of  this  technological  revolution. 

As  reported  in  the  fourteenth  edition  of  Ships  and  Aircraft  of  the 
U.  S.  Fleet  the  AV-8B  Harrier  II  VSTOL  attack  aircraft  "provides  twice 
the  payload  of  the  AV-8A  with  up  to  9200  lbs  of  external  stores."  The 
primary  difference  between  the  AV-8A  and  the  AV-8B  is  the  extensive 
use  of  composite  materials  in  the  wing  and  fuselage.  This  is  just  one 
example  of  the  use  of  composite  materials  in  existing  military  applica- 
tions. New  applications  include  the  VS-22  Osprey  tilt  rotor  aircraft, 
whose  wing  and  fuselage  will  be  entirely  made  of  graphite  epoxy  com- 
posite material.  Emerging  "stealth"  technology  and  its  application  to 
aircraft  is  a  further  new  application  of  composites. 

Composites  offer  designers  many  attractive  features.  These 
include  very  high  strength-to-weight  ratios,  high  stiffness-to-weight 
ratios,  wear  resistance,  controllable  heat  expansion  coefficients,  cor- 
rosion resistance,  fatigue  resistance,  low  price,  and  ease  of  manufac- 
turing and  repairability.  The  field  is  expanding  rapidly  and  research  is 


bringing  to  the  designer  stronger,  tougher,  and  stiffer  materials  every 
day. 

A  hurdle  to  be  overcome  for  composites,  however,  is  the  general 
lack  of  understanding  of  the  physics  and  mechanisms  at  work  in  the 
composite  which  give  these  materials  their  unique  properties.  This 
lack  of  understanding  has  restrained  designers'  ability  to  accurately 
predict  composite  strength.  To  compensate  for  this,  large  factors  of 
safety  are  being  used  in  the  design  of  composite  structures,  which  in 
turn  adds  weight  and  size  and  causes  other  related  problems  that 
reduce  design  efficiency.  As  the  application  of  composite  materials  has 
expanded  to  larger  and  larger  components,  such  as  entire  fuselages, 
rocket  motor  cases,  and  gun  turrets,  this  problems  has  become 
magnified. 

Composite  materials  are  two  or  more  dissimilar  materials  that 
when  properly  combined  form  a  new  material  whose  properties 
exceed  either  of  the  individual  constituents.  For  the  designer,  com- 
posite materials  provide  an  ability  not  only  to  optimize  his  design  but 
also  to  optimize  the  material  to  be  used  for  the  design.  Material  prop- 
erties of  conventional  materials  that  can  be  improved  upon  by  using 
composites  include  stiffness,  corrosion  resistance,  wear  resistance, 
strength-to-weight  ratio,  and  radar  wave  absorbability  and  reflectivity. 
Composites  in  use  today  include  fiberglass,  WEST  (Wood  Epoxy 
Saturation  Technique),  graphite-epoxy,  Aramid-epoxy,  and  Boron- 
epoxy.  Composite  materials  also  include  graphite -aluminum  and  hybrid 
composites  of  more  than  two  constituents. 


Commonly  accepted  types  of  composites  are  the  fibrous  compos- 
ites, which  consist  of  fibers  in  a  matrix;  laminated  composites,  which 
consist  of  layers  of  various  materials;  and  particulate  composites,  which 
are  composed  of  particles  in  a  matrix.  The  composite  is  fabricated  by 
combing  the  reinforcement  material  into  a  generally  ductile  matrix 
such  as  epoxy.  This  lay-up  is  called  a  lamina.  Several  layers  of  lamina 
may  be  combined  in  prescribed  geometrical  orientations  to  create  a 
laminate.  The  designer  specifies  his  materials  and  then  designs  the 
laminate  to  meet  his  requirements. 

Until  the  early  1980s  material  properties  of  composites  were 
predicted  by  a  simple  principle  called  the  "rule  of  mixtures."  This 
approximation  assumed  that  each  constituent  of  the  composite  con- 
tributed to  the  overall  composite  property  in  direct  proportion  to  its 
property  weighted  by  its  volume  percentage  in  the  composite.  This 
linear  weighting  method  was  satisfactory  for  composite  physical  prop- 
erties, which  are  averaged  combinations  of  the  constituent  properties 
such  as  density  and  stiffness. 

However,  the  "rule  of  mixture"  is  very  inaccurate  for  those  physi- 
cal properties  which  are  governed  by  extreme  values  such  as  perme- 
ation and  strength.  For  such  properties,  considerable  testing  of  a 
specific  composite  is  needed  to  achieve  representative  material  prop- 
erty values  suitable  for  engineering  design  and  reliability  assurance. 
Not  only  are  large  amounts  of  data  required  for  each  specific 
composite  but  the  designer's  material  selections  become  limited  to 
only  those  composites  which  have  been  tested.  This  characterization 


methodology    could    not    keep    up    with    the    rapid    pace    of    the 
introduction  of  new  materials. 

Clearly,  a  characterization  method  was  needed  to  be  able  to 
predict  composite  material  properties,  particularly  strength,  if  the  full 
advantage  of  composites  was  to  be  obtained.  Much  work  has  been 
completed  in  the  area  by  researchers  such  as  Phoenix,  Rosen,  Harlow, 
and  Wu.  The  most  current  papers  on  reliability  and  strength  predic- 
tions of  composites  may  be  found  in  References  1,  2,  3,  and  4.  The 
most  significant  part  of  their  efforts  is  that  they  have  identified  that 
the  strength  of  composites  is  a  probabilistic  manifestation  of  the  sta- 
tistical strength  of  the  fiber  strength.  The  statistics  of  the  fiber 
strength  govern  the  number  of  failure  sites  in  a  composite.  The  proba- 
bilistic spatial  distribution  of  the  failure  sites  together  with  the  failure 
mechanism  govern  the  strength  of  a  specific  composite. 

Different  probabilistic  models  have  evolved  to  capture  the  essence 
of  the  failure  mechanisms.  The  simplest  model  is  the  "weakest  link 
model,"  which  envisions  the  entire  composite  as  a  long  thin  fiber;  the 
entire  chain  (the  composite)  fails  when  the  weakest  link  of  the  fiber 
fails.  This  model  does  not  account  for  the  redundancy  provided  by  the 
matrix  and  as  a  result  grossly  underestimates  the  strength  of  the  com- 
posite. The  "equal  load-sharing  model"  accounts  for  the  redundancy 
provided  by  the  matrix  through  definition  of  an  effective  length,  an 
isolation  parameter  for  the  initial  breaks.  However,  this  model  does 
not  account  for  the  stress  concentration  around  the  breaks  and,  as  a 
result,    over-predicts    the    composite    strength.    Finally,    the    "local 


load-sharing  model"  accounts  for  both  the  matrix  redundancy  and  the 
local  stress  concentration.  These  refinements  lead  to  the  most 
realistic  predictions  to  date. 

There  are  two  phenomena  relating  to  the  statistical  strength  of 
composites.  The  first  is  that  the  larger  a  composite  becomes,  the 
weaker  it  gets.  The  second  phenomenon  is  that  the  larger  a  compos- 
ite, the  less  strength  scatter  among  such  large  structures.  This  is  due 
to  the  the  lower  weak  tail  of  the  strength  distributions  of  the  con- 
stituent fibers.  This  effect  has  been  demonstrated  in  the  laboratory 
and  in  the  field. 

An  extension  of  this  theory  then  becomes  that  if  the  fibers  are 
preloaded  to  break  the  weak  sections  of  the  fiber  prior  to  their  lay-up 
into  the  matrix,  then  the  lower  tail  of  the  composite  strength 
distribution  is  reduced  and,  in  fact,  becomes  bounded  by  a  value.  This 
process  is  called  prestressing.  For  a  prescribed  value,  the  strength  and 
reliability  of  the  composite  may  then  be  theoretically  calculated  for  a 
given  fiber  and  the  magnitude  of  preload  based  on  the  modified  fiber 
distribution.  The  variables  to  be  determined  to  prestress  a  composite 
include  the  magnitude  of  the  preload,  the  gauge  length  to  which  the 
load  is  applied,  and  the  modified  fiber  parameters. 

The  ability  to  accurately  recover  the  fiber  strength  parameters  for 
small  (less  than  2.5  cm)  gauge  lengths  is  very  important  if  the  pre- 
stress effect  is  to  be  optimized.  The  bundle  testing  done  in  the  labora- 
tory yields  data  that  includes  not  only  the  fiber  data  but  also  machine 
compliance,    slack,    and    friction.    For    short    gauge    lengths,    the 


contributions  of  slack  and  friction  are  small  but  compliance  becomes 
relatively  large.  The  portion  of  displacement  contributed  to  by  the 
experimental  system  compliance  in  a  bundle  test  needs  to  be  removed 
from  the  data  if  the  true  fiber  parameters  are  to  be  obtained. 

This  study  entailed  the  mathematical  simulation  modeling  of 
strength  prediction  for  a  bundle  of  fibers  with  no  matrix.  This  simula- 
tion was  run  on  an  IBM  Personal  Computer  using  Microsoft  Fortran 
4.01  for  source  code  and  Lotus  1-2-3  for  graphing.  When  the  simula- 
tion program  had  been  validated,  a  zero  gauge  length  bundle  test  was 
completed  to  measure  the  INSTRON  fiber  bundle  testing  machine 
compliance.  A  set  of  procedures  was  written  to  optimize  a  curve  fit  of 
the  compliance  data  to  provide  an  explicit  equation  which  models  the 
INSTRON  system  compliance. 

When  this  procedure  had  been  validated,  bundle  tests  of  varying 
lengths  were  tested  at  the  Mechanics  of  Materials  for  Composites  Lab- 
oratory at  the  Naval  Postgraduate  School.  The  displacement  con- 
tributed to  the  test  data  by  the  INSTRON  system  compliance  was 
removed  from  the  total  displacement.  Under  another  thesis  effort,  this 
rectified  data  was  analyzed  in  order  to  extract  the  fiber  statistical 
strength  distribution  parameters. 

The  purpose  of  this  investigation  was  to  validate  the  experimental 
procedures  and  data  analysis  methods  with  the  use  of  computer  simu- 
lated data.  Data  generated  from  Monte  Carlo  simulation  using  known 
parameters  is  then  utilized  to  assist  in  experimental  design  and  the 
associated    prediction    theory.     Experimental    design    and    data 
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interpretation  methodology  provide  a  proven  foundation  toward  iden- 
tification of  an  optimum  gauge  length  for  prestressing  in  order  to 
guarantee  the  reliability  of  composite  structures. 


II.     BACKGROUND 

The  tensile  failure  process  of  composite  materials  is  a  complex 
sequential  combination  of  many  different  processes.  These  processes 
are  controlled  by  the  statistical  strength  of  the  fiber  and  the  relative 
spatial  clustering  of  the  broken  fiber  sites  within  a  composite.  Given  a 
probabilistic  model  which  is  capable  of  characterizing  the  sequential 
failure  mechanisms,  a  formula  can  be  derived  to  predict  the 
probabilistic  strength  of  the  composite  in  terms  of  the  statistical 
strength  of  the  fiber.  The  parameters  for  such  a  probabilistic  models 
can  then  be  expressed  in  terms  of  other  common  measures  of  central 
tendencies  (e.g.,  mean  composite  strength)  and  dispersions  (e.g., 
strength  scatter). 

A.      COMPOSITE  STRENGTH  PREDICTION 

The  prediction  of  the  probabilistic  strength  of  a  composite  as  a 
function  of  fiber  statistics  was  first  introduced  by  B.  W.  Rosen  in 
Reference  5,  in  the  early  sixties.  He  refined  the  simple  weakest  link 
model  (which  is  essentially  a  series  model)  by  accounting  for  the 
redundancy  introduced  by  the  presence  of  a  matrix  binder.  The  com- 
posite geometry  he  examined  is  shown  below: 

•  a  single  layer  of  fibers 

•  evenly  distributed  throughout  the  matrix 

•  no  fibers  are  in  contact  with  other  fibers 

•  the  load  is  purely  axial  and  in  tension. 


8 


The  strength  of  a  single  lamina  of  composite  material  is  con- 
trolled by  four  factors: 

•  The  strength  and  modulus  distributions  of  the  fiber. 

•  The  strength  and  modulus  distributions  of  the  matrix. 

•  The  interaction  between  the  fiber  and  the  matrix  with  respect  to 
load  sharing  and  stress  concentration. 

•  The  location  distribution  of  relative  strengths  of  the  fibers  with 
respect  to  the  location  distribution  of  relative  strengths  of  the 
matrix. 

When  a  composite  is  pulled  in  tension,  because  the  fiber  generally 
has  a  much  higher  modulus  as  compared  to  the  matrix  (e.g.,  graphite 
108  Psi  versus  epoxy  104  psi),  the  fibers  will  carry  most  of  the  load. 
The  load-carrying  capability  of  the  matrix  can  be  neglected. 

As  the  load  is  increased,  each  fiber  will  carry  an  equal  amount  of 
load  that  is  equal  to  the  applied  load  divided  by  the  number  of  fibers. 

Pf  =  Pt/n  (1) 

The  load  on  each  fiber  will  be  equal,  but  because  the  diameter  of 
the  fibers  is  not  constant,  the  stress  carried  at  each  point  in  the  fiber 
will  be  different.  The  stress  will  vary  inversely  to  the  diameter  of  the 
fiber.  Some  of  the  variations  are  regions  of  gradual  thinning,  others  are 
indentations  or  nicks.  These  areas  of  reduced  thickness  or  indenta- 
tions or  nicks  are  called  "flaws."  These  flaws  reduce  the  load-carrying 
capability  of  the  fiber  by  reducing  its  thickness  and  introducing  stress 
concentration  areas  in  the  fiber. 


A  great  deal  of  study  has  been  done  in  identifying  and  quantifying 
the  effect  of  these  flaws  on  a  microscopic  scale.  The  work  of  Phoenix 
and  Harlow  proposed  that  on  a  macroscopic  scale,  if  the  flaws  are  dis- 
tributed randomly  with  respect  to  location  on  the  fibers  and  normally 
with  respect  to  the  severity  of  the  flaw,  then  the  strength  of  a  fiber 
can  be  modeled  as  a  long  chain  of  segments.  The  strength  of  the  seg- 
ments may  be  modeled  using  the  Weibull  distribution.  The  model 
selection  is  expanded  in  Appendix  A.  The  Weibull  distribution  is  an 
appropriate  model  for  the  probability  of  failure  of  a  chain  with  many 
links. 

F(x)  =  1  -  exp  I-(x/px)<x]  (2) 

Where      x      =  is  the  load  applied  to  the  fiber  (the  chain) 

a     =  is  the  shape  parameter  characterizing  the  variability 

|3X    =  is  the  scale  parameter  characterizing  central  location 
of  the  failure  strengths. 

The  Weibull  distribution  for  strength  predicts  the  "size  effect." 
That  is,  if  the  parameters  are  known  for  one  given  length,  the 
parameters  for  any  other  lengths  of  the  same  material  can  be  com- 
puted. This  is  also  expanded  in  Appendix  A. 

As  the  load  increases  on  the  material,  the  stress  on  each  segment 
of  the  fiber  will  increase  as  well.  When  the  load  reaches  a  point  where 
the  resulting  stress  is  equal  to  the  ultimate  strength  of  a  segment,  that 
segment  of  fiber  will  break.  This  break  will  release  a  certain  amount  of 
strain  energy  into  the  composite.  It  will  also  cause  an  immediate  local 
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load  concentration  that  will  be  felt  by  the  adjacent  and  surrounding 
fibers. 

The  matrix,  via  its  shear  stress  carrying  capability,  will  attempt  to 
absorb  the  released  strain  energy  and  distribute  the  local  load  concen- 
tration to  the  surrounding  fibers.  The  actual  values  of  the  stress  con- 
centration factors  caused  by  the  fiber  break  are  not  known,  but  shear 
leg  analysis  yields  reasonable  approximations. 

Once  one  fiber  segment  breaks,  many  events  can  occur.  The  first 
is  that  the  matrix  is  strong  enough  to  absorb  the  released  strain 
energy  and  the  surrounding  fiber  elements  are  capable  of  carrying  the 
local  load  concentration.  If  this  occurs,  the  material  is  stable  and  the 
load  on  each  fiber  will  increase  to 

Pf  =Pt/(l  -  n)  (3) 

Where  n  represents  the  number  of  broken  fibers. 

Another  possibility  is  that  the  matrix  can  not  absorb  the  released 
strain  energy.  In  this  case  the  matrix  will  split  as  it  fails  in  shear.  The 
material  will  continue  to  tear  until  its  radius  is  such  that  the  stress 
concentration  decreases  to  a  point  below  the  ultimate  strength  of  a 
section  of  the  matrix.  The  fibers  in  this  region  will  then  be  subjected 
to  additional  stress  as  the  matrix  pulls  away.  If  the  segments  of  the 
matrix  are  strong  enough  to  carry  the  additional  load,  the  material  will 
then  again  become  stable. 

A  further  possibility  is  that  the  adjacent  fiber  segment  may  not  be 
able  to  withstand  the  stress  caused  by  the  first  fiber  failure.  In  this 
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case,  that  fiber  segment  will  break  and  cause  additional  strain  energy 
to  be  released.  This  reaction  is  called  auto  clustering.  The  process  will 
continue  until  it  is  interrupted  by  a  fiber  segment  that  is  strong 
enough  to  carry  the  stress  or  the  entire  piece  fails  catastrophically. 

Some  observations  of  the  description  above:  The  composite  may 
fail  due  to  either  the  local  strength  of  the  matrix  or  the  local  strength 
of  the  fiber  or  a  combination  of  both.  The  likely  initiation  of  the  failure 
process  will  occur  with  a  fiber  break  because  the  strength  scatter  or 
variance  of  the  fiber  is  typically  greater  than  that  of  the  matrix  giving 
rise  to  nucleation  sites  at  already  high  stress  regions.  The  chain  of 
events  that  occur  after  that  first  break  is  dependent  on  the  strength 
statistics  of  the  adjoining  matrix  and  fibers.  If  there  are  many  strong 
segments  surrounding  a  weak  segment  that  fails,  then  additional 
material  will  not  fail.  If,  however,  this  is  not  the  case,  then  further 
damage  occurs  which  may  lead  to  final  unstable  catastrophic  failure. 

From  the  previously  described  failure  process,  it  can  be  observed 
that  the  strength  of  the  composite  is  dependent  on  the  nucleating 
weak  segments  and  the  probability  of  there  being  another  weak  seg- 
ment nearby.  Therefore,  the  shape  factor  of  the  fiber  strength  distri- 
bution is  the  critical  element  that  controls  composite  strength  and 
reliability.  Because  of  the  mechanism  of  failure,  as  a  composite 
becomes  larger,  the  number  of  the  population  of  fiber  segments 
increases  and  the  likelihood  of  the  existence  of  very  weak  segments 
becomes  larger.  This  causes  the  strength  of  the  entire  composite  to 
decrease. 


12 


Another  observation  is  that  because  the  strength  of  the  composite 
is  controlled  by  the  weak  fiber  segments,  the  strength  may  not  be 
solved  for  explicitly  because  the  actual  value  of  the  weak  segments  is 
probabilistic.  The  strength  of  adjoining  segments  is  probabilistic,  as 
well. 

Prestressing  solves  two  critical  problems.  First,  it  sets  the  lower 
limit  of  the  strength  of  the  segments  by  eliminating  the  weak  seg- 
ments by  causing  them  to  break  prior  to  the  application  of  matrix. 
Second,  prestressing  causes  the  weak  segments  to  be  broken  without 
local  load  stress  concentrations  affecting  the  adjacent  fibers  and  with- 
out the  localized  dynamic  strain  energy  release  because  the  weak  seg- 
ments are  broken  before  the  matrix  is  solidified.  This  results  in  auto 
clustering  being  minimized. 

The  critical  implementation  problem  is  then  to  find  the  "segment 
length"  to  ensure  that  the  majority  of  the  weak  segments  are  elimi- 
nated and  then  determine  a  prestress  level.  The  prestress  level  may 
be  optimized  to  a  level  where  the  advantages  of  breaking  the  weak 
segments  is  balanced  by  the  introduction  of  clusters  to  the  matrix. 

B.     PRESTRESSING 

G.  J.  Mills  and  associates,  under  an  United  States  Air  Force-spon- 
sored project,  conducted  prestress  tests  on  Boron  fibers  in  the  early 
1970s.  In  his  reports,  documented  in  References  5,  6,  and  7,  he 
detailed  how  his  work  demonstrated  the  feasibility  of  the  concept  and 
provided  data  from  early  experiments  that  supported  his  expectations. 
His    paper    did    not    attempt    to    model    or    predict    the    effect    of 
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prestressing.  He  did  note  the  effect  of  varying  prestress  gauge  length 
and  that  as  the  gauge  length  decreased,  the  strength  of  the  fiber 
increased,  but  he  did  not  attempt  to  quantify  an  optimum  length.  He 
described  the  strength  increase  of  the  composite  purely  in  terms  of 
removal  of  weak  segment  and  reducing  the  tail-end  scatter  of  the 
fiber.  He  did  not  attempt  to  predict  the  performance  of  the  composite 
once  the  matrix  had  been  set.  He  did  report  that  the  strength  did 
increase. 

In  1986,  Lt.  David  Bell  and  Professor  E.  Wu  of  the  Naval  Postgrad- 
uate School  conducted  prestress  tests  on  AS-4  graphite/epoxy.  The 
results  are  published  in  Reference  9.  The  report's  experimental  data 
confirmed  the  phenomenon  of  strength  improvement  after  pre- 
stressing. This  was  in  concurrence  with  Mills'  laboratory  data.  Both 
papers  accounted  for  the  improvement  by  theorizing  that  the  process 
of  prestressing  improved  the  strength  of  a  composite  because  it 
reduced  the  strength  scatter  of  the  composite.  The  work  done  by  Bell 
and  Wu  was  only  done  with  one  gauge  length,  10",  and  even  with  the 
problems  encountered  in  handling  the  samples,  the  material  did  show 
a  significant  reduction  in  scatter  and  therefore  an  increase  in  the 
usable  strength  reliability. 

Implementation  of  these  findings  to  a  practical  application 
requires  a  model  that  can  adequately  predict  the  effect  of  prestressing 
as  a  function  of  prestress  level,  gauge  length,  and  fiber  strength 
parameters. 
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C.     TESTING  METHODS 

There  are  two  methods  in  which  fiber  strength  statistics  can  be 
measured.  One  method  is  to  prepare  and  perform  a  breaking  strength 
test  on  one  individual  fibers.  The  second  method  is  to  form  many 
fibers  into  a  bundle  of  parallel  fibers  and  then  perform  a  controlled 
strain  rate  failure  test. 

The  first  method  is  very  time  consuming  and  error  prone.  The 
fibers  are  very  thin  (in  the  order  of  10  |im)  and  handling  requires 
delicate  treatment.  Laboratory  conditions  must  be  carefully  monitored 
as  the  small  fibers  are  susceptible  to  damage  caused  by  temperature, 
wind,  abrasion,  and  other  disturbances.  In  addition,  each  test  only 
provides  one  single  point  of  data.  If  failure  probability  data  is  required 
down  to  .001  failure  rate,  then  at  least  10,000  samples  should  be 
tested. 

Obviously,  the  second  method  is  the  preferred  method.  The  bun- 
dle test,  performed  in  our  case  on  bundles  of  1,000  filaments,  will 
yield  1,000  data  points  with  one  test.  The  handling  problems,  while 
still  present,  are  on  a  much  more  manageable  level.  The  data  analysis 
of  the  load  vs.  displacement  data  will  yield  the  strength  distribution 
parameters  for  the  fiber  that  is  equal  to  1,000  individual  filament 
tests. 

It  is  important  then  that  we  measure  the  true  load  and  displace- 
ment for  just  the  fibers.  The  displacement  that  the  testing  machine 
measures  is  the  sum  of  the  fiber  displacement,  machine  compliance, 
electronic  noise,  slack  in  the  individual  fibers,  and  friction  between 
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the  fibers.  Figure  1  demonstrates  that  test  data  is  the  combination  of 
all  of  these  factors. 


Test  Data 


Friction 


Figure  1 
Typical  Bundle  Test  Data  Graph 

The  compliance  is  due  to  the  mechanical  construction  of  the 
testing  machine.  The  machine  operates  with  gears,  clutches,  and 
motors.  Each  of  these  has  an  initial  start-up  slip  or  backlash.  The 
machine  also  has  grips  which  hold  the  samples  and  the  bases  which 
will  elongate  under  load.  They  are  very  stiff  and  the  elongation  is  small 
when  compared  to  the  bundles,  except  when  the  bundle  gauge  length 
is  short.  This  is  illustrated  in  Figure  2,  where  the  displacement  is 
shown  as  percent  strain.  The  compliance  for  the  long  sample  has  a 
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greater  slope  than  the  short  sample  because  the  slope  of  the  compli- 
ance in  displacement  multiplied  is  by  the  gauge  length  for  percent 
strain. 


compliance 


compliance 


bundle  test 


bundle  test 


%  Strain        3  %  Strain     3 

GL  =   Long  0  25  cm)  GL  =   Short  (<  2.5  cm) 

Figure  2 

Effect  of  Compliance  on  Long-Gauge  Length 
vs.  Short-Gauge  Length 

The  compliance  of  the  machine  can  be  measured  if  the  there  is  no 
contribution  to  the  displacement  from  the  fibers.  This  is  achieved 
using  a  zero  gauge  length  sample.  A  sample  must  be  utilized  to  include 
the  sample  grip  mechanism  in  the  compliance.  A  curve-fitting 
approximation  is  then  applied  to  the  data  to  create  a  continuous  func- 
tion of  d  in  terms  of  P.  This  function  may  then  be  subtracted  from  the 
test  data  to  eliminate  the  effect  of  compliance.  This  procedure  and 
testing  is  expanded  in  Appendix  C. 
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The  effect  of  noise  is  sensed  in  the  load  cell.  The  sensitivity  of  the 
load  cell  is  1  out  of  10,000.  This  is  due  to  the  limit  of  the  analog-to- 
digital  converter,  which  can  only  manage  four  bits  of  data.  The  amount 
of  noise  in  the  machine  sensed  by  the  load  cell  is  very  small  when 
compared  to  the  sensitivity  of  the  load  cell.  As  the  load  increase  and 
the  output  changes  in  order  of  magnitude,  the  noise  is  not  even 
noticeable.  The  noise  in  the  machine  is  considered  to  be  random 
electrical  white  noise  due  to  interference  in  the  data  transmission 
lines.  Observation  of  the  data  indicates  that  the  noise  is  not  very 
significant. 

A  simulation  program  was  written  during  this  study  that  allowed 
the  user  to  input  random  noise  up  to  any  level.  This  will  allow  simula- 
tions to  be  run  and  investigate  the  effect  of  noise  on  the  fiber  strength 
parameters  recovered  from  the  simulated  test. 

The  slack  is  due  to  the  limitation  of  the  sample  preparation  that 
all  the  fibers  do  not  have  exactly  the  same  gauge  length.  This  problem 
is  found  to  be  greater  with  the  longer  samples  as  compared  to  the 
short  samples.  The  slack  can  not  be  eliminated  by  current  sample 
preparation  techniques,  but  extreme  care  by  the  technician  in 
preparing  the  sample  can  minimize  the  effect. 

Lt.  Joseph  Schmidt  is  examining  the  effects  of  slack  on  the  calcu- 
lated fiber  parameters  and  expects  to  publish  the  results  of  his  study 
in  September  1988,  from  the  Naval  Postgraduate  School,  Monterey. 

The  effect  of  friction  in  the  bundle  test  cannot  not  be  definitively 
characterized.  The  presence  of  friction  is  manifested  by  large  vertical 
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drops  in  the  load  displacement  curve,  perhaps  due  to  the  tangling  of 
broken  fibers  with  unbroken  fibers.  The  friction  also  will  cause  a  fibers 
to  break  more  than  once.  The  effects  of  friction  do  increase  with 
longer  gauge  lengths  and  are  reduced  with  shorter  gauge  lengths. 

With  these  factors  that  affect  the  data  of  the  testing  machine 
accounted  for,  the  true  fiber  strength  statistics  may  be  recovered  from 
the  bundle  tests.  These  statistics  may  then  by  used  to  optimize  the 
prestress  effect  by  calculating  the  effect  prestressing  will  have  in 
altering  the  fiber  strength  distribution. 
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III.  PRESTRESSING  MECHANISM 

The  prestressing  process  is  to  preload  the  composite  fibers  prior 
to  solidification  of  the  matrix  to  cause  the  weak  segments  of  the  fibers 
to  break.  The  load  that  was  carried  by  the  broken  filaments  is  then 
shared  equally  among  the  remaining  fibers.  Equal  load  sharing  is  in 
effect  because  the  pre-solidified  matrix  cannot  support  a  shear  load 
and  cause  stress  concentration  cells  to  form.  Upon  solidification  of  the 
matrix,  however,  local  load  sharing  is  in  effect,  thereby  increasing  the 
local  redundancy.  Because  the  weak  segments  of  the  fibers  were 
broken  prior  to  solidification  of  the  matrix,  the  local  load  concentra- 
tion stress  cells  are  not  as  severe  as  if  the  fibers  had  broken  while  in 
the  matrix.  The  reduction  in  the  severity  of  the  stress  concentration 
cells  provides  the  reliability  enhancement  of  the  composite  and  is  the 
motivation  for  prestressing. 

The  parameters  that  affect  the  prestress  effect  are  gauge  length, 
prestress  magnitude,  and  the  fiber  strength  distribution  statistics.  The 
optimum  gauge  length  to  prestress  the  fibers  would  be  the  ineffective 
length.  The  ineffective  length,  as  described  by  Rosen  in  Reference  5, 
is  a  function  of  both  the  tensile  strength  of  the  fiber  and  of  the  shear 
strength  of  the  matrix.  This  makes  the  ineffective  length  unique  for 
each  composite,  and  because  both  strengths  are  random  variables,  the 
ineffective  length  for  a  particular  location  in  the  composite  is  random 
as  well. 
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As  the  gauge  length  increases,  the  likelihood  of  breaking  all  of  the 
weak  segments  in  the  fiber  decreases.  This  is  because  using  the  equiv- 
alent load  sharing  model,  it  is  only  possible  to  cause  one  break  in  each 
fiber  over  the  gauge  length.  The  break  occurs  at  the  weakest  flaw, 
which  may  occur  before  reaching  the  desired  prestress  level.  As  illus- 
trated in  Figure  3,  this  allows  other  flaws  weaker  than  the  prestress 
level  to  remain  unbroken. 
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Figure  3 
Effect  of  Preloading  for  Various  Gauge  Lengths 
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When  the  prestressed  fiber  is  placed  into  a  matrix,  the  local  load- 
sharing  model  becomes  applicable.  With  this  model,  the  weak  flaws  in 
the  composite  will  break  under  load  and  generate  the  larger  stress 
concentration  cells  which  prestressing  is  designed  to  minimize.  The 
mathematical  modeling  of  solving  for  the  probability  density  function 
of  remaining  weak  flaws  given  gauge  lengths  greater  than  the  ineffec- 
tive length  and  prestress  levels  is  very  complex.  This  problem  is  being 
treated  in  a  separate  project. 

The  preload  level  has  two  effects  on  the  prestress  effect.  For 
higher  prestress  levels,  more  fiber  breaks  are  introduced  and  more 
stress  concentration  cells  will  be  created  in  the  the  composite  once 
the  matrix  has  solidified.  The  severity  of  the  concentration  cells  is  a 
function  of  the  ineffective  length.  The  shorter  the  ineffective  length  is, 
the  smaller  the  stress  concentration  cell  is  because  of  less  local  load 
sharing.  Second,  the  higher  the  preload  level,  the  more  difficult  it 
becomes  to  grip  or  maintain  the  desired  gauge  length.  The  fibers  by 
their  nature  have  a  small  diameter,  on  the  order  of  5  microns  for 
graphite,  and  are  very  susceptible  to  handling  damage. 

The  underlying  fiber  statistics  also  have  an  effect  on  the  prestress 
effect.  The  fiber  strength  variability  (the  Weibull  model  shape  parame- 
ter a),  has  the  greatest  influence  on  the  prestress  effect.  As  seen  in 
Figure  4  ,  if  the  shape  parameter  is  very  high  (a  >  20  ),  i.e.,  the  fiber 
has  low  variability,  then  the  prestress  effect  will  not  be  very  helpful. 
This  is  because  the  range  of  the  distribution  of  the  fiber  strengths  is 
very  small  and  the  weak  tail  that  prestressing  eliminates  or  reduces  is 
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inherently  not  present.  Fibers  of  this  strength  uniformity,  however, 
are  not  available  at  this  time.  If  the  shape  parameter  is  small,  then  the 
scatter  is  much  greater.  With  this  scatter,  the  weak  tail  may  be  altered 
or  removed  by  prestressing. 
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Figure  4 
Strength  Distribution  for  Fibers  (Weibull  Model) 

The  quantification  of  the  prestressing  effect  can  be  written  as  a 
function. 


PS  =  f  (O,  ¥,  a,  p) 


(4) 


where     O  =  Gauge  Length 

¥  =  Preload  Level 

a  =  Fiber  shape  parameter,  Weibull 

p  =  Fiber  scale  parameter,  Weibull 
The  functional  interrelation  between  these  functions  is  presently 
not  known  and  appropriate  probabilistic  modeling  will  be  required. 


23 


One  approach  is  to  optimize  each  variable  separately  where  possible. 
For  example,  the  Weibull  strength  parameters  characterize  the  fiber 
strength  distribution.  The  prestress  effect  of  the  preload  level  and  the 
gauge  length  will  alter  these  parameters.  The  optimization  of  the 
effect  on  the  alteration  of  these  parameters  on  the  reliability  of  the 
composite  can  be  investigated  as  one  means  of  optimizing  the  pre- 
stress effect. 

The  probability  of  failure  (F*)  vs  normalized  load  (x*/(3)  graph  for 
the  fibers  (F)  and  composite  (C)  with  the  Failure  Probability  Density 
Function  superimposed  is  shown  in  Figure  5.  This  figure  demonstrates 
how  the  PDF  can  be  transposed  to  a  linear  curve.  The  linear  curve  is 
much  easier  to  use  and  because  it  reveals  in  great  clarity  what  is  hap- 
pening at  the  weak  tail.  This  graph  will  be  used  to  show  the  effect  of 
prestressing  with  respect  to  failure  probability. 


Fiber  Statistics, f  +  Failure  Model 
=>  Composite  Failure  Probability 


For  Composite  Service  Risk,  Fs 
Fiber  Statisitcs  Must  be  Assured 
for  x  >  xc 


Normalized  strength 

Figure  5 
Probability  of  Failure  vs.  Normalized  Strength 
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The  probability  of  failure  for  the  composite  (C)  is  calculated  from 
the  fiber  parameters.  The  composite  reliability  will  diverge  from  the 
fiber  reliability  at  Xc  because  of  the  local  load-sharing  stress-concen- 
tration cells.  These  cells  cause  the  composite  fibers  in  the  matrix  to 
break  earlier  because  the  fibers  adjoining  a  broken  fiber  have  to  carry 
the  load  carried  by  broken  fiber  by  themselves,  not  spread  equally 
amongst  all  of  the  remaining  fibers.  Until  the  load  Xc,  the  two  curves 
do  not  diverge  and  the  reliability  of  the  composite  is  equal  to  the  reli- 
ability of  the  fibers. 

If  the  fibers  have  been  damaged  so  that  their  distribution  curve  is 
modified  to  appear  like  that  shown  in  Figure  6,  the  result  on  the  com- 
posite strength  can  be  seen.  As  the  composite  graph  swings  upward, 
its  reliability  is  decreasing.  The  composite  is  stronger  than  the  fiber 
due  to  the  load  sharing  but  its  reliability  for  loads  greater  than  Xc  has 
diminished.  The  composite  reliability  for  loads  below  Xc  has  not 
changed  at  all. 

If  the  fibers'  strength  distribution  could  be  improved,  the  failure 
probability  curve  would  be  shifted,  as  shown  in  Figure  7.  The  impact 
on  the  composite  strength  is  that  the  curve  swings  downward.  This 
downward  swing  reflects  greatly  increased  reliability.  In  fact,  the 
composite  will  not  fail  theoretically  at  all  for  any  load  less  than  Xc. 

A  method  to  achieve  the  fiber  distribution  in  Figure  7  would  be  to 
proof  test  all  of  the  composite  fibers  up  to  Xc  and  then  remove  any 
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Figure  6 
Probability  of  Failure  vs.  Normalized  Strength  (Damaged  Fibers) 


Normalized  strength 


Figure  7 
Probability  of  Failure  vs.  Normalized  Strength  (Improved  Fibers) 


26 


broken  fibers  from  the  bundle  prior  to  adding  the  matrix.  The  effect  of 
this  process  on  the  failure  probability  curve  is  shown  in  Figure  8  as 
curve  R.  The  cost  of  this  operation,  considering  the  extreme  thinness 
of  the  fibers,  the  susceptibility  to  damage  due  to  handling,  and  the 
extreme  amount  of  fibers  that  make  up  the  composite,  make  this 
process  in  reality  impossible. 


X 


r 

R 


NRL 


-  Proof  test  Fibers 

"  Original  distribution  f 

-  Remove  broken  fibers 
Resultant  distribution 

-  Not  Remove  broken  fibers 
Long  Gauge  Length 
Resultant  distribution 


NRS  -  Not  Remove  broken  fibers 
Short  Gauge  Length 
Resultant  distribution 


X 


Figure  8 
Effects  on  Failure  Probability  Curves  Due  to  Proof  Testing 

Instead,  what  may  be  done  is  to  preload  the  fibers  to  Xc  on  a  gauge 
length  equal  to  the  ineffective  length  of  the  fiber.  Because  of  the  han- 
dling difficulties,  the  broken  fibers  will  not  be  removed.  The  effect  of 
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this  process  on  the  failure  probability  curve  is  shown  in  Figure  8  as 
curve  NRS. 

The  NRS  curve  has  swung  downward  as  compared  to  the  original 
fiber  curve,  but  it  has  not  swung  as  far  down  as  if  the  broken  fibers  had 
been  removed.  The  reason  it  does  not  swing  as  far  is  the  local  load- 
sharing  stress-concentration  cells  that  are  created  when  the  matrix  is 
added. 

If  the  fibers  are  preloaded  to  Xc  but  on  a  gauge  length  that  is 
greater  than  the  ineffective  length,  the  effect  on  the  failure  probability 
curve  is  shown  on  Figure  8  as  curve  NRL.  This  curve  swings  downward 
from  the  original  fiber  curve  even  less  because  of  the  probability  that 
some  fibers  remain  with  segments  that  are  weaker  than  the  preload. 

As  previously  stated,  the  strength  of  the  composite  may  be  calcu- 
lated from  the  strength  of  the  fibers.  If  the  the  fiber  strength  distribu- 
tions are  altered  as  shown  in  Figure  8,  the  composite  strength 
distribution  will  be  altered  in  the  same  direction  as  well.  Using  this,  a 
designer  of  a  composite  may  then  use  a  failure  probability  curve  to 
determine  whether  his  design  meets  his  reliability  specifications.  If  it 
does  not,  he  may  then  decide  to  improve  the  composite  by  prestress- 
ing  the  fibers  to  a  critical  load  to  alter  the  failure  probability  to  come 
within  the  requirements. 

The  movement  of  the  failure  probability  curve  can  be  predicted  by 
probability  modeling  based  on  the  preload  level,  the  gauge  length,  the 
ineffective  length,  and  the  original  fiber  strength  parameters.  The 
results  of  such  an  analysis  (by  Wu  and  Harlow)  on  the  fiber  swing  for  a 
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fiber  preloaded  at  the  level  0.3  *  p  ,  for  an  a  of  5.0  and  prestressed  at 
multiples  of  its  ineffective  length,  are  shown  in  Figure  9.  The  graph 
demonstrates  that  even  if  the  prestress  gauge  length  is  very  large 
compared  to  the  ineffective  length,  significant  improvement  of  fiber 
strength  distribution  can  be  achieved. 

The  optimization  of  the  prestress  effect  is  to  identify  the  ineffec- 
tive length  or  to  find  the  shortest  prestress  gauge  length  that  can  be 
applied. 
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Figure  9 

Probability  of  Failure  vs.  Normalized 
Strength  Prestress  Effect  Simulation 
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IV.   EXPERIMENTATION 

The  experimental  work  was  directed  in  two  areas.  The  first  was  to 
measure  the  machine  compliance  of  the  INSTRON  4206  (  a  universal 
materials  testing  machine)  in  order  to  characterize  the  testing 
machine  system  compliance  denoted  by  the  function  8C  (P).  Once  this 
was  completed,  the  second  portion  of  the  experimentation  was  to 
perform  failure  tests  on  AS-4  Graphite  bundle  samples  of  various 
lengths  in  order  to  determine  the  strength  characteristics  and  to 
observe  the  effects  of  gauge  length  on  the  test  results  with  respect  to 
compliance,  slack,  noise,  and  friction.  Bundle  samples  were  tested 
with  no  oil  and  then  with  oil  in  an  attempt  to  reduce  the  effects  of 
friction,  particularly  on  the  longer  samples  of  50  and  25  cm.  Small- 
gauge  lengths  of  2.5  and  1  cm  were  tested  with  particular  interest 
with  respect  to  the  examination  of  their  suitability  to  preloading. 

A.      COMPLIANCE   TESTING 

The  compliance  testing  was  accomplished  using  a  sample  illus- 
trated in  Figure  10.  The  copper  tabs  were  prepared  and  washed  with  a 
dilute  acid  to  enhance  the  performance  of  the  adhesive.  The  graphite 
bundle  was  glued  to  the  lower  set  of  copper  tabs  with  acetylene  adhe- 
sive. The  upper  tabs  were  then  set  so  that  the  distance  (d)  between 
the  the  upper  and  lower  tabs  was  as  near  zero  as  possible.  The 
graphite  bundle  was  then  glued  into  place  with  the  same  adhesive.  The 
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sample   was   prepared   by  Composites   Laboratory  Technician,   Jim 
Nageotte. 


Adhesive 

Graphite 
Bundle 

Adhesive 


Copper  Tabs 


d  =  0 


Copper  Tabs 


Figure  10 
Compliance  Test  Sample 

The  tensile  test  was  performed  on  the  INSTRON  4206  testing 
machine  and  associated  INSTRON  software,  version  4.01  .A  plot  of  the 
load  vs.  displacement  is  shown  as  Figure  11.  The  data  and  test  sum- 
mary forms  generated  by  the  software  are  enclosed  in  Appendix  E. 

The  compliance  data  was  then  fit  to  the  function  shown  as  Equa- 
tion 4.  This  was  done  in  accordance  with  the  procedure  for  curve  fit- 
ting detailed  in  Appendix  C.  The  curve  fitting  was  accomplished  with 
the  use  of  LOTUS  1-2-3  spreadsheet  program  and  an  IBM/AT.  The 
resulting  coefficients  for  the  curve  fit  of  the  compliance  are  listed  in 
Table  I. 
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Figure  11 
Compliance  Test  Load  vs.  Displacement  Graph 
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INSTRON  4206  Compliance  Curve  Fit  Function 


For  0  <  5  <  Dc  1 
For  Del  <  5 


5  (P)  =Ai(P)n  +  A2 
5  (P)  =  As(P)  +  A4 


(4) 


TABLE  1 
INSTRON  4206  COMPLIANCE  CURVE  FIT  COEFFICIENTS 


Coefficient 

Value 

Units 

Ai 

.0512 

mm 

A2 

-.0217 

mm 

n 

.1975 

— 

A3 

.0064 

mm /kg 

A4 

.0242 

mm 

B.     BUNDLE  TESTING 

The  bundle  testing  was  accomplished  with  the  use  of  samples  as 
illustrated  in  Figure  12.  The  copper  tabs  were  prepared  and  washed 
with  a  dilute  acid  to  enhance  the  performance  of  the  adhesive.  The 
graphite  bundle  was  glued  to  the  lower  set  of  copper  tabs  with  acety- 
lene adhesive.  The  upper  tabs  were  then  set  so  that  the  distance  (d) 
between  the  the  upper  and  lower  tabs  was  the  target  gauge  length. 
The  graphite  bundle  was  held  straight  with  a  2  kg  weight  while  the 
upper  tabs  were  glued  in  place.  The  samples  were  carefully  stored  in 
the  Composites  Laboratory  until  use.  As  in  the  case  of  the  compliance 
test  sample,  all  of  the  samples  tested  were  prepared  by  Laboratory 
Technician,  Mr.  Jim  Nageotte. 
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Figure  12 
Graphite  Bundle  Sample 

The  bundle  tests  were  performed  on  AS-4  Graphite  bundles  of 
3,000  filaments.  The  tests  were  run  in  accordance  with  the  procedure 
outlined  in  Appendix  D.  The  displacement  due  to  compliance  was  sub- 
tracted from  the  total  displacement  and  a  summary  of  the  the  results 
is  listed  in  Table  2.  The  data  and  test  summary  forms  generated  by  the 
software  are  enclosed  in  Appendix  E. 
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TABLE  2 
BUNDLE  TESTING  RESULTS  (COMPLIANCE  REMOVED) 

Material:      AS-4  Graphite  (3000  filaments) 


Test  # 

Dry/ 

Gauge 

Modulus 

a 

P 

Comments 

Oil 

Length 
(kg/mm) 

(mm/mm 

) 

090901 

Dry 

5.0 

cm 

1068 

3.214 

.0135 

090902 

Dry 

.5 

cm 

* 

090903 

Dry 

.5 

cm 

* 

090904 

Dry 

.5 

cm 

* 

090905 

Dry 

.5 

cm 

* 

090906 

Dry 

.5 

cm 

* 

090907 

Dry 

5.0 

cm 

* 

090908 

Dry 

5.0 

cm 

912 

4.465 

.0145 

090910 

Dry 

2.5 

cm 

* 

080901 

Oil 

70.0 

cm 

Test  Validatioi 

100901 

Dry 

2.5 

cm 

915 

3.000 

.0157 

100902 

Oil 

2.5 

cm 

* 

100903 

Oil 

2.5 

cm 

* 

100904 

Dry 

50.0 

cm 

905 

3.035 

.0077 

100905 

Oil 

50.0 

cm 

959 

3.660 

.0085 

100906 

Oil 

50.0 

cm 

Not  Analyzed 

100907 

Oil 

25.0 

cm 

932 

3.964 

.0101 

100908 

Dry 

25.0 

cm 

917 

3.180 

.0097 

*  Adhesive  Failure 

Examination  of  the  load  vs.  percentage  strain  plots  revealed  the 
fact  that  for  the  shorter  gauge  lengths,  2.5  and  1  cm,  the  copper  tabs 
failed  to  grip  the  graphite  bundle  uniformly.  This  caused  some  of  the 
fibers  to  slip  out  of  the  tabs  and  the  data  from  these  tests  were  of  no 
validity.  The  tests  that  had  this  phenomena  are  listed  as  adhesive  fail- 
ure in  the  comments  section  of  Table  2.  Figure  13  shows  a  represen- 
tative example  of  this  phenomenon  as  the  load  remains  constant  while 
the  percent  strain  increases  beyond  the  limit  of  the  experiment. 
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Figure  13 
Load  vs.  Percent  Strain  for  Sample  100903  (Adhesive  Failure) 
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An  examination  of  Table  2  shows  that  seven  of  the  18  samples 
resulted  in  satisfactory  tests.  The  seven  successful  tests  results  have 
been  summarized  into  Table  3. 

TABLE  3 

BUNDLE  TESTING  RESULTS  (COMPLIANCE  REMOVED) 

SUCCESSFUL  TESTS 

Material:      AS-4  Graphite  (3000  filaments) 


Test  # 

Dry/ 

Gauge 

Modulus 

a 

P 

Oil 

Length 
(kg/mm) 

(mm/mm) 

100901 

Dry 

2.5  cm 

915 

3.000 

.0157 

090901 

Dry 

5.0  cm 

1068 

3.214 

.0135 

090908 

Dry 

5.0  cm 

912 

4.465 

.0145* 

100908 

Dry 

25.0  cm 

917 

3.180 

.0097 

100907 

Oil 

25.0  cm 

932 

3.964 

.0101 

100904 

Dry 

50.0  cm 

905 

3.035 

.0077 

100905 

Oil 

50.0  cm 

959 

3.660 

.0085 

*  Data  not  plotted 

Figures  14  through  32  are  the  graphs  of  load  vs.  percent  strain  of 
the  test  data,  of  the  compliance,  and  of  the  test  data  with  the  fiber 
removed. 


38 


rs 

0 

0 

J 


AS-4  BUNDLE  TEST   TEST  DATA 

SAMPLE  10O901    GAJGE  LENGTH  =  25  MM 


0.5  1.2  1.E 

Percent  Strain  X 


Figure  14 
Load  vs.  Percent  Strain  Sample  100901  Test  Data 
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AS-4-  BUNDLE  TEST    COMPLIANCE 

SAfcPlI  100901    GAXJE  LENCTH  =  25  MM 
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0.5 
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Percent  Strain  X 
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Figure  15 
Load  vs.  Percent  Strain  Sample  100901  Compliance 
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AS -4  BUNDLE  TEST   COMPLIANCE  REMOVED 
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Figure  16 
Load  vs.  Percent  Strain  Sample  100901  Compliance  Removed 
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AS-4  BUNDLE  TEST   TEST  DATA 

SAkPLE  090301    OUGE  LDMGrTH  =  50  MM 


0.4  0.S  1.2  1.B 

Percent  Strain  % 


2.4 


Figure  17 
Load  vs.  Percent  Strain  Sample  090901  Test  Data 
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AS-4  BUNDLE  TEST   COMPLIANCE 

SWuPLE:  090901    {ACE  LENGTH  =  50  UW 


Percent  Strain  X 

Figure  18 
Load  vs.  Percent  Strain  Sample  090901  Compliance 
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AS-4-  BUNDLE  TEST   COMPLIANCE  REMOVED 
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SAJyPLE  Q3M01    OUGE  LENGTH  -  50  MM 


Percent  Strain  X 

Figure  19 
Load  vs.  Percent  Strain  Sample  090901  Compliance  Removed 
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AS -4  BUNDLE  TEST   TEST  DATA 

SAJuPLE  100308  GAUGE  LENGTH  =  250  Ml 
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Figure  20 
Load  vs.  Percent  Strain  Sample  100908  Test  Data 
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AS-4-  BUNDLE  TEST   COMPLIANCE 

SAMPLE  100906  GAUGE  LENGTH  =  250  Hi 


1.2 

Percent  Strain  X 


Figure  21 
Load  vs.  Percent  Strain  Sample  100908  Compliance 
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AS-4  BUNDLE  TEST   COMPLIANCE  REMOVED 
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Figure  22 
Load  vs.  Percent  Strain  Sample  100908  Compliance  Removed 
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AS-4  BUNDLE  TEST  (OIL)  TEST  DATA 

9WPH  10M07  GAUGE  LENGTH  =  250  MM 


Percent  Strain  % 


Figure  23 
Load  vs.  Percent  Strain  Sample  100907  Test  Data 
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AS-4  BUNDLE  TEST  (OIL)  COMPLIANCE 

SAMPLE  1GM07  GAUGE  LENGTH  =  250  MM 


Percent  Strain  X 


Figure  24 
Load  vs.  Percent  Strain  Sample  100907  Compliance 
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AS-4  BUNDLE  TEST  (OIL)  COMPLNC  REMOVED 
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Figure  25 
Load  vs.  Percent  Strain  Sample  100907  Compliance  Removed 
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AS-4  BUNDLE  TEST   TEST  DATA 

SUFI!  100904  GAUGE  LENGTH  =  SCO  Muf 
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Figure  26 
Load  vs.  Percent  Strain  Sample  100904  Test  Data 
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AS-4  BUNDLE  TEST   COMPLIANCE 

SAIuPLE  100904  GAUGE  LENGTH  =  500  W 
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Figure  27 
Load  vs.  Percent  Strain  Sample  100904  Compliance 
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AS-4  BUNDLE  TEST   COMPLIANCE  REMOVED 
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Figure  28 
Load  vs.  Percent  Strain  Sample  100904  Compliance  Removed 
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AS-4  BUNDLE  TEST  (OIL)  TEST  DATA 

SNPLE  10Q905  GAUGE  LENGTH  =  500  kM 


1.2  1.6 

Percent  Strain  X 


Figure  29 
Load  vs.  Percent  Strain  Sample  100905  Test  Data 
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AS-4  BUNDLE  TEST  (OIL)  COMPLIANCE 

SAAPLE  10Q905  GAUGE  LENGTH  =  500  Ml 
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Figure  30 
Load  vs.  Percent  Strain  Sample  100905  Compliance 
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AS-4  BUNDLE  TEST  (OIL)  COMPLNC  REMOVED 


SAMPLE  100905  GAUGE  LENGTH  =  £00  MM 


0.4  O.S  1.2  1.6 

Percent  Strain  X 


2.4 


Figure  31 
Load  vs.  Percent  Strain  Sample  100905  Compliance  Removed 
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V.  CONCLUSIONS 

A.      COMPLIANCE   TESTING 

The  curve-fitting  procedure  outlined  in  Appendix  C  was  successful 
in  arriving  at  a  continuous  function  to  describe  the  load-induced  dis- 
placement of  the  load  train  for  the  INSTRON  testing  machine.  Figure 
32  shows  the  data  points  measured  during  the  zero  gauge  length  test. 
Figure  33  shows  the  curve  fit  function  using  the  parameters  listed  in 
Table  2.  Figure  34  then  shows  the  data  points  and  the  curve  fit 
superimposed  on  one  another.  The  greatest  deviation  in  load  is 
.0017  kg. 

Based  on  the  small  deviation  of  the  data  points  to  the  curve  fit 
function,  the  curve  fit  of  the  compliance  is  satisfactory  to  be  used  to 
subtract  the  displacement  due  to  compliance  from  bundle  test  data. 

The  curve  fit  solution  for  the  machine  compliance  is  satisfactory 
for  loads  greater  than  .1  kg.  The  curve  fit  does  not  do  a  good  job  in  the 
region  below  this  load.  This  is  due  to  the  inconsistent  zero  displace- 
ment length.  The  zero  point  for  displacement  is  very  difficult  to 
achieve  for  bundle  testing.  The  sensitivity  of  the  load  cell  is  insuffi- 
cient to  consistently  identify  an  actual  zero  load  level  on  the  bundle. 
The  effect  of  this  on  the  fiber  strength  distribution  parameters  is  that 
some  elements  of  the  lower  tail  are  distorted. 
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Figure  32 
0  Gauge  Length  INSTRON  Compliance  Test 
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Figure  33 
INSTRON  Compliance  Curve  Fit  Function 
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Figure  34 
0  Gauge  Length  INSTRON  Compliance  Test  and  Curve  Fit 
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The  compliance  testing  method  for  the  INSTRON  has  been 
demonstrated  to  be  feasible  and  to  provide  realistic  results.  Testing 
should  continue  to  determine  whether  the  compliance  is  repeatable  or 
whether  it  is  a  function  of  other  factors  such  as  temperature,  humid- 
ity, or  age  of  the  machine.  Full  understanding  of  compliance  is  impor- 
tant because  of  the  magnitude  of  compliance  effects  on  short  gauge 
lengths;  it  is  the  short  gauge  lengths  to  which  stressing  is  applied. 

B.     BUNDLE  TESTING 

The  bundle  tests  demonstrated  the  difficulty  in  obtaining  test  data 
for  short  gauge  lengths.  Table  2  illustrated  the  fact  that  no  successful 
1  cm  tests  were  completed  and  only  one  2.5  cm  test  was  completed. 
An  improved  technique  of  sample  preparation  to  prevent  the  slipping 
of  the  fibers  through  the  tabs  is  needed.  Possibly  an  adhesive  with  bet- 
ter wetting  qualities  could  be  tried. 

The  difficulty  in  testing  the  small  gauge  lengths  highlights  two 
important  effects  on  prestressing.  The  first  is  that  the  shorter  the 
prestress  gauge  length,  the  more  difficult  the  process  becomes  due  to 
the  difficulty  in  maintaining  gauge  length.  The  second  important  effect 
is  that  the  prestress  effect  on  fiber  parameter  improvement  will  have 
to  be  based  on  a  gauge  length  less  than  the  ineffective  length.  The 
optimization  process  will  now  include  another  variable,  the  minimum 
possible  gauge  length  due  to  machine  constraints.  As  the  minimum 
gauge  length,  increases  the  effectiveness  of  the  prestress  will 
decrease. 
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The  importance  of  accurately  removing  the  displacement  in  the 
test  data  for  short  gauge  lengths  was  demonstrated  by  the  bundle 
tests.  Figures  15,  18,  21,  24,  27,  and  30  show  that  as  a  function  of 
percent  strain,  the  shorter  the  gauge  length,  the  greater  the  percent- 
age of  test  displacement  is  due  to  machine  compliance.  Figure  15 
shows  that  the  percent  strain  displacement  due  to  compliance  for  the 
shortest  gauge  length  (2.5  cm)  was  greater  than  .6  percent,  or  nearly 
25  percent  of  the  total  test  displacement.  Figure  27  shows  that  the 
percent  strain  displacement  due  to  compliance  for  the  longest  gauge 
length  (50  cm)  was  less  than  .02  percent,  which  is  less  than  1  percent 
of  the  total  test  displacement. 

The  inaccuracy  incurred  in  not  eliminating  the  compliance  dis- 
placement for  long  samples  before  reducing  the  data  to  determine  the 
fiber  strength  parameters  is  minimal.  In  the  case  of  short  samples, 
however,  this  is  not  the  case.  If  the  compliance  displacement  is  not 
removed,  the  effect  on  the  strength  parameters  will  be  the  same  as  if 
the  gauge  length  was  actually  longer  than  what  was  tested.  This  results 
in  the  fiber  parameters  being  conservative  in  the  prediction  of  fiber 
strength.  While  this  is  safe,  it  defeats  the  purpose  of  prestressing. 

Figures  25  and  31  are  the  graphs  of  the  oiled  bundles.  Comparing 
these  graphs  to  Figures  22  and  28,  which  are  the  dry  or  not  oiled 
samples,  the  effects  of  friction  can  be  seen  in  the  shape  of  the  curves 
after  the  maximum  load  has  been  reached.  The  dry  sample  load  drops 
in  large  vertical  jumps.  These  large  jumps  correspond  to  many  fibers 
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breaking  at  once  instead  of  sequentially.  Figure  35  shows  the 
simulated  test  data  for  a  100  mm  sample.  The  curve  is  very  smooth 
and  does  not  show  the  sharp  vertical  drops  the  actual  test  data  shows. 

The  large  vertical  drops  in  the  actual  test  data  are  noticeably 
reduced,  but  not  eliminated,  in  the  oiled  samples.  The  oil,  a  soybean 
derivative,  reduces  the  friction  but  does  not  eliminate  all  of  the  effects. 
The  effect  of  friction  on  the  test  does  appear  to  dominate  the  tail  end 
or  the  final  30  percent  of  the  percent  strain.  The  effect  of  friction  on 
the  fiber  strength  parameters  is  difficult  to  predict.  One  effect  that  is 
clear  is  that  the  data  for  the  first  50  percent  to  60  percent  of  the 
breaks  in  the  long  samples  do  not  appear  to  be  affected  by  friction  as 
the  remaining  breaks.  In  addition,  Figures  16  and  19  indicate  that  the 
effect  of  friction  declines  as  the  gauge  length  becomes  shorter. 

The  need  to  eliminate  friction  between  the  fibers  in  a  bundle  test 
is  important  if  long-gauge  length  tests  are  to  be  accurately  completed. 
Different  oils  may  be  tried  as  lubricants  to  improve  upon  the  perfor- 
mance of  the  soybean  oil.  If  the  friction  cannot  be  eliminated,  it  must 
be  brought  down  to  a  minimal  level. 

The  effect  of  noise  in  the  bundle  tests  appears  to  be  minimal.  Fig- 
ures 36,  37,  and  38  are  graphs  of  simulated  bundle  tests  with  noise 
added  to  the  load  data.  Figure  4  is  a  graph  of  the  same  simulation  with 
no  noise  present. 

A  noise  level  of  1  gm  is  equivalent  to  .002  percent  for  the  50  kg 
load  cell  used  by  the  INSTRON.  The  difference  in  the  curves  from 
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Figure  35 
Bundle  Test  Simulation  (Gauge  Length— 10  cm) 
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Figure  36 
Bundle  Test  Simulation  With  ±1  gm  Noise 
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Figure  37 
Bundle  Test  Simulation  With  ±10  gm  Noise 
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Figure  38 
Bundle  Test  Simulation  With  ±100  gm  Noise 
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Figure  35  to  Figure  36  is  barely  noticeable.  A  noise  level  of  10  gms  is 
equivalent  to  .02  percent.  Figure  37  shows  some  effects  of  noise,  but 
they  are  minimal.  A  noise  level  of  100  gms  is  equivalent  to  .2  percent. 
This  is  an  excessive  amount  of  noise  and  the  effects  of  a  noise  level 
this  high  are  clearly  shown  by  Figure  8. 

The  INSTRON  manual  specifies  the  sensitivity,  or  noise  level,  of 
the  load  cell  to  be  .01  kg.  This  corresponds  to  roughly  the  noise  level 
of  Figure  37.  An  inspection  of  Figure  37  demonstrates  that  the  effect 
of  this  noise  level  is  barely  discernable  when  compared  to  Figure  35,  a 
simulation  with  no  noise. 

This  research  has  demonstrated  that  the  physical  process  of 
preloading  fibers  with  a  short  gauge  length  is  difficult  due  to  the 
physical  limitations  of  the  sample  gripping  mechanism.  More  work  is 
required  to  find  an  adequate  means  to  grip  or  restrain  the  fibers  for 
short  gauge  lengths. 

The  removal  of  the  displacement  due  to  compliance  in  the  test 
results  is  significant  for  gauge  lengths  shorter  than  2.5  cm  in  order  to 
determine  accurately  the  fiber  strength  distribution  parameters.  The 
effect  of  compliance  diminishes  with  longer  gauge  lengths.  The  fric- 
tion between  the  fibers  in  bundle  alters  the  shape  of  the  load  vs.  dis- 
placement curve.  Friction  has  more  of  an  effect  on  longer  samples 
than  on  short  samples.  Friction  may  be  reduced  by  using  oil  on  the 
samples,  but  it  does  not  entirely  eliminate  the  effect.  Noise  is  not  a 
significant  factor  in  altering  the  data  and  does  not  interfere  with  the 
strength  distribution  parameters. 
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VI.   RECOMMENDATIONS 

Further  study  into  the  magnitude  of  the  machine  compliance  for 
the  INSTRON  4206  testing  machine  should  be  initiated.  The  effects  of 
temperature,  the  gripping  devices,  and  repeatability  should  be  investi- 
gated. To  complete  this,  a  quicker  method  to  determine  the  curve  fit 
coefficients  needs  to  be  developed.  A  Fortran  program  that  can  accu- 
rately and  quickly  solve  the  curve-fitting  routine  in  Appendix  C  should 
be  written  and  validated.  A  series  of  compliance  tests  may  then  be 
completed  to  build  a  database. 

The  methods  for  sample  preparation  and  bundle  testing  need  to 
become  less  technician  dependent.  There  are  too  many  areas  where 
mishandling  of  the  samples  may  introduce  errors  into  the  data  by 
damaging  the  fibers.  The  reliance  on  the  laboratory  technician  to  build 
the  samples,  align  the  samples,  and  then  set  the  gauge  length  all  add 
variability  to  the  data  that  should  be  reduced  if  not  eliminated. 
Hydraulic  grips  can  eliminate  the  hand-  and  wrench-tightening  pro- 
cedure required  by  the  current  grips.  A  brace  or  stand  should  be 
designed  and  built  to  hold  the  sample  during  loading  into  the 
INSTRON.  This  will  eliminate  much  of  the  handling  the  sample  must 
now  endure. 

The  prestress  effects  of  preloading  the  fibers  to  a  load  greater 
than  the  critical  load  should  be  investigated.  This  is  a  probability 
problem  whose  solution  could  be  integrated  into  the  bundle  simulation 
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program.  This  program  could  then  be  used  to  predict  the  prestressing 
performance  and  then  validated  with  tests  performed  in  the  laboratory 
on  gauge  lengths  that  can  be  successfully  tested.  This  is  very  important 
because  the  poor  test  results  of  the  short  samples  indicate  that  a  short 
preload  gauge  length  (<  1  cm)  may  be  difficult  to  achieve. 

The  effects  of  prestressing  with  respect  to  long-term  reliability 
and  composite  fatigue  failure  should  be  investigated.  The  impact  of 
local  load  stress  cells  on  long-term  reliability  needs  to  be  addressed 
quickly.  The  urgency  of  this  aspect  of  prestressing  is  that  due  to  the 
nature  of  fatigue  failure,  data  taking  is  a  very  long  process. 
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APPENDIX  A 
MODEL   SELECTION 

Composite  fiber  materials  such  as  graphite,  aramid,  or  boron  are 
characteristically  long  and  thin  when  compared  to  the  dimensions  of 
the  composite.  For  example,  the  diameter  of  a  graphite  fiber  is  on  the 
order  of  5  microns,  while  the  fiber  length  can  extend  to  hundreds  of 
feet  for  a  1  square  foot  pressure  vessel.  Under  tension,  the  fibers  fail, 
or  break,  at  positions  where  flaws  or  imperfections  exist.  These  flaws 
cause  the  fiber  to  break  at  loads  much  below  the  theoretical  strength 
of  the  fiber,  which  is  the  chemical-bonding  strength.  The  type  of  the 
flaw  may  be  scratching,  chemical  erosion,  contamination,  or  one  of 
many  others.  The  source  of  the  flaw  may  have  come  in  production, 
handling,  or  service  use.  The  spatial  location,  density,  and  severity  of 
these  flaws  determine  the  strength  of  each  fiber. 

A  single  fiber  or  filament  may  be  visualized  as  a  chain  made  up  of 
many  links.  These  links,  or  segments,  each  have  a  unique  strength 
that  is  limited  by  the  severity  of  the  flaws  contained  in  that  segment. 
The  more  severe  the  flaw,  the  weaker  the  segment.  For  a  long  fiber, 
one  with  many  segments,  the  strength  of  that  fiber  is  only  the 
strength  of  the  most  severe  flaw,  or  the  weakest  link.  If  the  severity  of 
the  flaw  is  known  to  fall  within  a  certain  probability  distribution  and 
the  location  of  the  flaws  are  randomly  distributed  along  the  fiber,  then 
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a  "weak  link"  model  can  be  assumed  for  predicting  the  strength  of  a 
fiber. 

From  statistics,  the  Weibull  distribution  best  fits  the  "weak  link" 
model.  Phoenix  and  Wu,  in  Reference  2,  describe  the  Weibull  statisti- 
cal model  as  it  applies  to  composite  fiber  strength  prediction  in  detail. 
In  general,  the  fiber  is  partitioned  into  a  series  of  segments  (m).  The 
shortest  limiting  length  segment  is  equal  to  the  ineffective  length  or 
effective  load  transfer  length  of  the  fiber.  This  is  on  the  order  of  mm 
for  graphite /epoxy.  A  bundle  is  then  formed  by  joining  several,  (n), 
fibers  in  parallel.  For  example,  in  this  project,  AS-4  bundles  where  m 
equals  1  and  n  equals  3,000  were  tested  in  the  laboratory. 

Based  on  the  Weibull  weak  link  model,  the  failure  cumulative  den- 
sity function  (CDF)  of  a  fiber  can  be  explicitly  written  as  a  function  of 
the  two  Weibull  parameters  a  and  p.  The  Weibull  failure  CDF  is  written 
as  Equation  5. 

F  =  1  -  exp  {-(x/P)<*}  (5) 

The  two  parameters,  a  and  p,  correspond  to  the  distribution 
shape  and  scale,  a  adjusts  the  shape  of  the  distribution  by  altering  the 
skew  and  range  of  values.  As  shown  in  Figure  39,  a  high  shape  para- 
meter, or  a,  corresponds  to  a  very  narrow  distribution;  an  a  equal  to 
3.5  approximates  a  normal  distribution;  and  a  smaller  a  corresponds 
to  a  very  wide  distribution  with  large  tails.  When  a  is  large,  the  scale 
parameter  p  is  approximately  the  mean.  The  AS-4  that  was  tested  in 
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this  project  for  a  gauge  length  of  25  cm  has  an  accepted  a  equal  to  4 
and  a  (3  equal  to  16  gm/ (mm/mm). 


P  X 

Figure  39 
Effect  of  the  Shape  Parameter  on  the  Weibull  Distribution 

The  Weibull  distribution  is  not  only  a  physically  appropriate  model 
for  fiber  strength,  it  has  an  additional  advantage  in  that  its  failure 
cumulative  distribution  (CDF)  can  be  written  as  an  explicit  equation.  In 
addition,  the  CDF  may  be  linearized  with  respect  to  p.  This  results  in 
Equations  6  and  7.  These  two  equations  may  then  be  graphed  together 
to  form  a  graph  such  as  the  one  shown  as  Figure  40.  This  graph  of  the 
F*/x*  function  allows  direct  graphical  correlation  of  strength  to  failure 
probability. 
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Figure  40 
Weibull  Linearized  Failure  Prediction  Graph 

F*  =  ln(-ln(l  -  exp  {-(P/p)a))) 
x*  =  In  (P/P) 


(6) 
(7) 


where  P  is  trie  load  subjected  to  the  filament. 

To  predict  reliability,  which  is  one  minus  the  failure  probability 
for  a  fiber,  the  a  and  (3  must  be  known  for  that  fiber  to  a  sufficient 
degree  of  accuracy.  Fiber  bundle  failure  tests  are  performed  in  the  lab- 
oratory to  arrive  at  the  parameters.  A  difficulty  arrives  with  the  fact 
that  the  parameters  are  a  function  of  the  length  of  a  fiber.  This  corre- 
sponds to  the  fact  that  a  long  chain  is  weaker  than  a  short  chain 
because  the  longer  chain  has  a  higher  probability  of  having  a  weak 
segment. 
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The  Weibull  model  provides  a  means  to  normalize  the  parameters 
between  different  lengths  of  a  common  fiber  if  the  parameters  are 
known  for  any  one  length.  This  allows  for  the  laboratory  results  to  be 
applied  to  fibers  of  lengths  different  than  the  lengths  tested.  This  fea- 
ture of  the  Weibull  function  is  commonly  called  the  size  effect. 

If  the  strength  distribution  parameters  pi  and  al  are  known  for  a 
fiber  of  length  li  and  the  parameters  p2  and  oc2  are  desired  for  the 
same  fiber  but  of  a  different  length  I2  the  following  derivation  provides 
the  relationship: 

Fi  =  1  -  exp  {-(x/pi)ai)  (5) 

R=  1  -F 
(Reliability  =  1  -  Failure)  (8) 

Ri  =  exp  {-(x/Pi)ai}  (9) 

m  =  li/l2  (10) 

RT  =  riRi  R2 =  (Ri)m 

RT  =  [exp  {-(x/P!)ai}]m  (11) 

R2  =  exp  {-m(x/pi)ai}  (12) 

R2  =  exp  {-{x/p2)a2}  (9) 

Equating  Equations  9  and  12  then: 
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exp  {-m(x/pi)ai}  =  exp  {-(x/p2)a2} 
{-m(x/p!)ai}  =  {-(x/p2)a2} 
(x/dn)"1/"!  (3i)ai  =  (x/p2)a2  (13) 

Solving  Equation  13  results  in  Equations  14  and  15. 

(X2  =  ai  (14) 

p2  =  PiCm)"1/*!  (15) 

An  implication  of  the  size  effect  relations,  Equations  14  and  15,  is 
that  a  long  composite  made  up  of  many  segments  may  be  reduced  to  a 
composite  of  one  segment.  This  simplification  provides  for  easier 
mathematical  modeling  of  the  composite  as  a  large  composite  model 
shown  in  Figure  41  with  many  segments  may  be  reduced  to  the  much 
smaller  model  shown  in  Figure  42.  This  greatly  reduces  the  amount  of 
computational  effort  required  to  model  composite  strength  or  bundle 
strength  as  the  model  matrix  has  been  reduced  from  a  (m  X  n)  to  a  (1 
Xn). 
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Figure  41 
Composite  Model  Without  Size  Effect 
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Figure  42 
Composite  Model  With  Size  Effect 
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APPENDIX  B 
BUNDLE  FAILURE  SIMULATION 

The  purpose  of  the  bundle  failure  simulation  was  twofold.  The  first 
was  to  allow  the  user  to  readily  produce  simulated  bundle  test  data 
with  variable  fiber  strength  distribution  parameters.  The  resulting  data 
could  then  be  graphically  presented  into  load  versus  displacement  or 
load  verses  percent  strain  representations.  These  graphs  were  utilized 
to  gain  understanding  and  insight  into  the  effect  of  fiber  strength 
parameters  and  gauge  length  on  fiber  failure  process. 

The  second  purpose  of  the  bundle  failure  simulation  was  to 
provide  a  means  of  verification  for  data  interpretation  and  reduction 
software  being  developed  for  use  in  conjunction  with  the  INSTRON 
Testing  machine  used  in  the  Composites  Laboratory,  Naval 
Postgraduate  School.  The  simulation  could  be  used  to  ensure  that  the 
software  returned  the  parameters  from  the  simulated  test  data  that 
were  used  as  inputs  to  the  bundle  simulation.  For  this  project,  the 
compliance  removal  procedures  and  software  described  in  Appendices 
C  and  D  were  validated  with  this  program.  Lt.  Joseph  Schmidt,  in  a 
separate  project,  used  this  program  to  validate  software  that  reduced 
bundle  test  data  into  the  fiber  strength  probability  distribution 
parameters. 

The  bundle  failure  simulation  program  "T2"  takes  fiber  strength 
distribution   parameters   a  and  P  for  the  length  of  the  bundle,  the 
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number  of  filaments,  the  length  of  the  bundle,  the  fiber  modulus,  and 
testing  machine  settings  of  cross-head  speed  and  sampling  rate,  and 
then  produces  two  sets  of  load  vs.  displacement  data.  The  first  set  is 
the  loaded  and  displacement  determined  for  each  individual  filament 
break.  This  is  referred  to  as  the  analog  data.  The  second  set,  which 
simulates  fiber  testing  machine  data,  is  the  continuous  load  and 
incremental  displacement  data,  which  is  referred  to  as  discrete  data. 
The  program  also  has  the  added  feature  of  inputting  a  noise  range.  The 
maximum  number  of  filaments  is  1,000. 

The  program  solves  for  the  strength  of  each  filament  based  on  the 
Weibull  distribution  and  the  size  effect.  An  explanation  of  the  Weibull 
function  and  of  the  size  effect  is  detailed  in  Appendix  A. 

A.     ANALOG  DATA  SIMULATION 

The  analog  data  was  determined  by  using  a  random  number 
generator  to  to  generate  a  random  strength  of  each  fiber,  determine 
the  displacement  of  each  break,  and  then  compute  the  load  of  the 
bundle  at  each  break. 

Equation  16  was  used  to  determine  the  strength  of  each  fiber. 

Pf  =  ln(-ln(l  -  exp  {-(x/p)"}))  (16) 

where    Pf  =  failure  load  of  the  the  fiber 
x  =  random  number  0  <  x  <  1 
Equation  17,  which  is  a  manipulation  of  Hooke's  law,  is  used  to 
determine  the  displacement  of  each  fiber  at  the  failure  load. 
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df=Pf*l/E  (17) 

where    df  =  Displacement  of  the  Fiber 

1  =  Length  of  the  Fiber  (Gauge  Length) 

E  =  Strain  Modulus 
The  displacement  at  each  break  was  stored  in  an  array  and  then 
sorted  in  ascending  order.  The  load  of  the  bundle  at  each  break  was 
then  determined  by  using  Equation  18.  This  is  also  a  manipulation  of 
Hooke's  Law.  At  each  break,  because  the  simulation  is  for  a  controlled 
displacement,  the  load  in  the  bundle  drops  by  a  factor  of  1/n.  This 
data  point  is  solved  for  by  Equation  19. 

(Pbh  =  ((dfh*  E*  (n-i  +1))/1  (18) 

(Pb)i+  =  ((dfh*  E*   (n-i))/l  (19) 

where    i  =  1,2,3,  .  .  .  n  (Number  of  Breaks) 

B.     DISCRETE  DATA  SIMULATION 

The  discrete  data  was  found  by  utilizing  the  analog  data.  The 
analog  data  was  transformed  into  a  continuous  function  by  connecting 
the  data  points.  The  displacement  between  data  points  is  found  by 
Equation  20. 

Deld  =  CHS/SR  (20) 
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where  CHS  =  Cross-head  Speed  (mm/min) 
SR  =  Sampling  Rate  (pts/sec) 
The  discrete  data  displacement  is  in  increments  of  Deld  while  the 
discrete  data  load  is  the  intercept  of  the  analog  continuous  data  for 
that  appropriate  displacement.  Equation  21  is  used  to  compute  the 
discrete  displacement  points.  Discrete  data  points  are  solved  for  until 
the  load  reaches  zero. 

d  =  Deld  *  I  1  =  0,1,2 Load  =  0  (21) 

C.  NOISE  SIMULATION 

The  user  has  the  option  of  simulating  the  effects  of  noise  in  the 
load  transducer  of  a  bundle  test.  If  desired,  the  noise  is  added  to  the 
discrete  load  data  by  Equation  22. 

Pnoise  =  P  -  «Ph  -  Pl)/2)  +  X  *  (Ph  -  Pi)  (22) 

where    P  =  Original  Load 

Ph  =  High  Noise  Range 
Pi  =  Low  Noise  Range 

D.  PROGRAM  LISTING 

The  program  was  written  in  Microsoft  Fortran  4.01.  It  is  designed 
to  be  operated  from  any  IBM /Compatible  computer  with  two  floppy 
drives  or  hard  disk  and  one  floppy  drive.  The  executable  program  is 
run  from  the  A:  drive  or  the  hard  disk  and  the  data  is  written  to  the  B: 
drive.  To  execute  the  program  type:  T5  <enter>. 
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SUBROUTINE  LISTING: 


INIT 

INPUT 

RAND 

SLOAD 

SORT 

ANLDAT 
DISDAT 

NOISE 

OUTPUT 


Initializes  arrays  to  zero 

Prompts  user  to  input  required  value 

Random  number  generator  (user  must  input  seed) 

Solves  for  the  breaking  stress  of  each  filament 

Sorts  the  filaments  by  breaking  stress  and  then 
solves  for  breaking  displacement 

Creates  the  analog  load  and  displacement  file 

Solves  for  the  continuous  and  creates  the  discrete 
load  and  displacement  file 

Allows  user  to  input  random  noise  to  a  specified 
threshold  into  the  discrete  load  data 

Writes  the   analog  and  discrete   data  to  floppy 
disks 


VARIABLE  LISTING: 

DATA(0: 1000,7)         Workspace  array 


ADATA(0:2000,2) 

DDATA(0:2000,2) 

ALPHA 

BETA 

E 

L 

N 

CHS 


Analog  data  output  file;  Displacement,  Load,  listed 
sequentially  by  first  break  to  last 

Discrete   data  output  file;   Displacement,   Load, 
listed  as  a  function  of  sampling  time 

Weibull  shape  parameter  for  the  given  length 

Weibull  scale  parameter  for  the  given  length 

Modulus  of  fiber 

Length  of  the  Fiber 

Number  of  filaments  in  the  bundle  (max.  =  1,000) 

Simulated  testing  machine  cross-head  speed 
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SR 
DELD 
POINTS 
SEED 

SEED1 

SEED2 

BREAKS 

LIMIT 

PH 

PL 


Simulated  machine  sampling  rate 

Distance  between  data  points  for  discrete  data 

Number  of  discrete  data  points 

Random    number   generator   seed    for    strength 
computations 

Random  number  generator  seed  for  discrete    data 
load  during  break  generator 

Random    number    generator    seed    for    noise 
generator 

Counter  that  is  used  to  find  slope  of  analog  data 
for  up  to  n- 1  breaks 

Maximum  strength  of  fiber  i  used  in  solving  for 
the  maximum  displacement  for  n-1  breaks 

Upper  noise  range 

Lower  noise  range 


OUTPUT  FILES 

ADATA 
DDATA 


Analog  data;  displacement,  load 
Discrete  data;  index,  displacement,  load 
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PROGRAM  T2 

PROGRAM  T2  WILL  SIMULATE  A  BUNDLE  TEST  FOR  N  FIBERS   IN  PARALEL 
FOR   GIVEN  FIBER   STRENGTH   PARAMETERS   AND   LENGTH   OF   THE   BUNDLE. 
THE   PROGRAM  WILL  OUTPUT  ANALOG  DATA  THAT   IS   A   LISTING   OF  THE 

*  DISPLACEMENT  AND   LOAD   OF  EACH  FIBER   BREAK.         THE   PROGRAM  WILL 

*  OUTPUT  DISCRETE   DATA  THAT   SIMULATES   MACHINE  TESTING  DATA  FOR  THE      * 

*  BUNDLE.       THE   OUTPUT   IS   BY   INDEX,    DISPLACEMENT  AND   LOAD.       THE 
NUMBER  OF  DATA   POINTS    IS   A  FUNCTION   OF  THE   SIMULATED   CROSSHEAD 
SPEED  AND  THE   SAMPLING  RATE. 

*  Vc 

*  THE   OUTPUT   IS  WRITTEN  TO  DRIVE   B:      ADATA.      AND  DDATA. 

•k  it 

*  DECLARATION  STATEMENTS 

REAL  DATA(0; 1000,7),    DDATA( 0; 5000 ,2) , ALPHA, BETA, L, CHS ,PH,PL, 
1E,DELD,SR,ADATA(0;  2000,2) 
INTEGER  N,I,J,K,POINTS,Q 
1   CONTINUE 

*  PRINT  BANNER  AND   START  PROGRAM 

PRINT*, 'COMPOSITE   BUNDLE   TEST   SIMULATION   PROGRAM1 

CALL   INIT( DATA, DDATA, ADATA) 

CALL   INPUT( ALPHA, BETA, L, CHS. E,N,SR) 

CALL  RAND ( DATA, N) 

CALL  SLOADC DATA, N, ALPHA, BETA) 

CALL   SORT(DATA,N,E,L) 

CALL  ANLDAT(DATA,N) 

CALL  DISDAT( DDATA, DATA, N,E,L, CHS, SR, POINTS) 

CALL  NOISE (DDATA, POINTS) 

CALL  OUTPUTC DATA , DDATA , N , E , L , ALPHA , BETA , POINTS , ADATA) 

*  ALLOW  USER  TO  RERUN   PROGRAM 

PRINT* 

PRINT*, 'RUN   COMPLETE,    DO   YOU  WISH  TO  RUN  ANOTHER   SET  OF  DATA* 

PRINT"--,'    IF   YES   TYPE    1,    ELSE   TYPE    0   TO   EXIT   PROGRAM' 

READ*,    Q 

IF(Q   .  EQ.  1)    GOTO    1 

*  PRINT   BANNER  AND  END   PROGRAM 

PRINT* 

PRINT*,'    END  OF  RUN,  HAVE  A  NICE  DAY' 

END 

SUBROUTINE   INIT  * 

SUBROUTINE   INIT   INITIALIZES  THE  ARRAYS  TO   ZERO  * 

y  .  y-  y  *  y^.  y-  y-  y-  y -  y-  y-  y-  y-  y-  y  «  y«  »t-  »t*  y  *  y.  y.  jl  .'-  „'..  j-  y-  - 1-  »>.. .*..  «•-  ju  -»*  j„  «',  »t„  j^  j„  y-  «».  ,«.  _f,  »f,  jl  y.  -.'-  -'.  y  .  «'-  ju  -'*  -*-  ju  y-  j-  y-  -•-  -j-  y-  J-  -*-  y.  y*  »'-  y*  y-  y  -  y-  y  -  y-  y-  y-  y- 

*  DECLARATION   STATEMENTS 

SUBROUTINE   INIT( DATA , DDATA , ADATA ) 

REAL  DATA( 0; 1000 , 7 ) ,DDATA( 0; 5000 , 2) , ADATA( 0; 2000 ,2) 

INTEGER   I, J 

*  PRINT   BANNER   AND   START   SUBROUTINE 
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PRINT* 

PRINT*,'*****  RUNNING  INITIALIZATION  SUBROUTINE  *****' 
PRINTS- 
DO  10  I  =  0,1000 
DO  9  J  =  1,7 
DATA(I,J)  =  0. 
9  CONTINUE 
10  CONTINUE 

DO  20  I  =  0,  5000 
DO  19  J  =  1,2 
DDATA(I.J)  =  0. 

19  CONTINUE 

20  CONTINUE 

DO  30  I  =  0,  2000 
DO  29  J  =  1,2 
ADATA(I,J)  =  0 

29  CONTINUE 

30  CONTINUE 
END 

SUBROUTINE  INPUT 

*  * 

SUBROUTINE  INPUT  PROMPTS  THE  USER  TO  INPUT  THE  REQUIRED  * 

DATA  TO  PERFORM  THE  SIMULATION.   INPUTS  ARE: 
v.-  ,'. 

ALPHA,  BETA,  N,  L,  E  -   FIBER  PARAMETERS 

CHS,   SR  -  TESTING  MACHINE  PARAMETERS  * 

*  Vr 


*  DECLARATION  STATEMENTS 

SUBROUTINE  INPUT( ALPHA , BETA , L , CHS , E , N , SR) 
REAL  ALPHA, BETA, CHS, L,SR,E 
INTEGER  N 

*  PRINT  BANNER  AND  START  SUBROUTINE 

PRINT* 

PRINT*,'*****  RUNNING  INPUT  SUBROUTINE  ****** 

PRINT* 

PRINT*, 'INPUT  ALPHA  AND  BETA  FIBER  PARAMETERS' 

READ*,  ALPHA,  BETA 

PRINT*, 'INPUT  NUMBER  OF  FILAMENTS  IN  BUNDLE  (MAX  =  1000)' 

READ*,  N 

PRINT*,' INPUT  MEAN  BUNDLE  LENGTH  (CM)' 

READ*,  L 

L  =  L*10 

PRINT*, 'INPUT  LOAD  MODULUS  OF  FIBER  (GM/MM/MM)' 

READ*,  E 

PRINT*, 'INPUT  SIMULATION  CROSSHEAD  SPEED  (MM/MIN)* 

READ*,  CHS 

PRINT*, 'INPUT  SIMULATION  SAMPLING  RATE  (PTS/SEC)1 

READ*,  SR 

PRINT* 

PRINT*, 'INPUT  COMPLETE' 

PRINT* 
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END 

*********************************************************************** 

*  SUBROUTINE  RAND 

*  * 

SUBROUTINE  RAND  GENERATES  A  RANDOM  NUMBER  FROM  0  TO  1  FOR  EACH  * 

FILAMENT  IN  THE  BUNDLE.   THE  RANDOM  NUMBER  GENERATOR  IS  STARTED  * 

FROM  THE  INITIAL  INPUT  SEED.   THE  FUNCTION  RAND  IS  USED  TO  * 

*  ACTUALLY  FIND  THE  RANDOM  NUMBER.  * 

*  * 
********************************************* ************************** 

SUBROUTINE  RAND(DATA,N) 

*  DECLARATION  STATEMENTS 

REAL  DATA(0; 1000, 7), RNG 
INTEGER  N,  SEED,  I 

*  PRINT  BANNER  AND  START  SUBROUTINE 

PRINT* 

PRINT*,'*****  RUNNING  RANDOM  NUMBER  GENERATOR  SUBROUTINE  ****** 
PRINT* 
PRINT* 

PRINT*,' INPUT  ANY  ODD  INTEGER  (STRENGTH  RNG  SEED)' 
READ*,  SEED 
PRINT* 

DO  100  I  =  1,  N 
DATA(I,1)  =  I 
DATA (I, 3)  =  RNG (SEED) 
100  CONTINUE 

END 
*********************************************************************** 

SUBROUTINE  SLOAD 

*  * 

SUBROUTINE  SLOAD  WILL  GENERATE  A  BREAKING  STRENGTH  FOR  EACH 
FILAMENT  AND  FROM  THE  MODULUS  FIND  THE  CORESPONDING  DISPLACEMENT.  * 

*********************************************************************** 

SUBROUTINE  SLOAD(DATA,  N, ALPHA, BETA) 

*  DECLARATION  STATEMENTS 

REAL  DATA(0; 1000, 7), ALPHA, BETA 
INTEGER  N 

*  PRINT  BANNER  AND  START  SUBROUTINE 

PRINT* 

PRINT*,'*****  RUNNING  SIMULATION  LOAD  SUBROUTINE  ****** 

PRINT* 

DO  100  I  =  1,  N 

*  BREAKING  STRENGTH  EQUATION  BASED  ON  WEIBULL  ALPHA  AND  BETA  PARAMETERS 

DATA( I ,4)  =  EXP( (LOG( -LOG( 1-DATA( 1,3)) )+ALPHA*LOG( BETA) ) /ALPHA) 
100  CONTINUE 
END 
*********************************************************************** 

SUBROUTINE  SORT 

*  * 

SUBROUTINE  SORT  USES  A  BUBBLE  SORT  METHOD  TO  SORT  THE  FILAMENTS 
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FROM  LOWEST  TO  HIGHEST  BY  BREAKING  POINT  DISPLACEMENT.  * 

}V  V- ^Wr  Vr  Vr  Vc  Vr  >WoWnW."V  Vr  V w  c  V- VoV -,V  iV  Vr  Vc  V^ 

SUBROUTINE  SORT(DATA,N,E ,L) 

*  DECLARATION  STATEMENTS 

REAL  DATA(0; 1000,7) ,T4,T5,E,L 
INTEGER  N, PAIRS ,T1 
LOGICAL  DONE 

*  PRINT  BANNER  AND  START  SUBROUTINE 

PRINT* 

PRINT*,1*****  RUNNING  SORT  SUBROUTINE  ****** 

PRINT* 

PAIRS=  N-l 

DONE  =  .FALSE. 
20  IF(.NOT.  DONE)  THEN 

DONE  =  .TRUE. 

DO  30  I  =  1, PAIRS 

IF(DATA(I,4)  .  GT.  DATA( 1  +  1,4))  THEN 

T4  =  DATA(I,4) 

Tl  =  DATA(I,1) 

DATA(I,4)  =  DATA( 1+1,4) 

DATA(I,1)  =  DATA( 1+1,1) 

DATA( 1+1,4)  =  T4 

DATA( 1+1,1)  =  Tl 

DONE  =  .FALSE. 

END  IF 
30  CONTINUE 

PAIRS  =  PAIRS  -  1 

GOTO  20 

ENDIF 

PRINT* 

PRINT*,'***  SORT  COMPLETE,  NEW  INDEX  BEING  ADDED  TO  DATA***' 

PRINT* 

DO  40  I  =  1,  N 

DATA(I,2)  =  I 
40  CONTINUE 

PRINT* 

PRINT*,'***  SOLVING  FOR  DISPLACEMENT  BASED  ON  BREAK  LOADS  ***' 

PRINT* 

DO  50  I  =  1,N 

DATA(I,5)  =  (DATA(I,4)*L)/(E) 
50  CONTINUE 

END 

ycVryrVryfy?yr7Vy?ynVyoV-WryoWc-WwrV«V^ 

*  SUBROUTINE  ANLDAT  * 

*  SUBROUTINE  ANLDAT  WRITES  THE  DISPLACEMENT  AND  LOAD  OF  EACH  BREAK   * 

*  INTO  ARRAY  AD ATA 

VoVyryoWoWrynWryfycyoWnVynWcy.-yoVVryry-yoVyf^ 

SUBROUTINE  ANLDAT( DATA , N ) 

*  DECLARATION  STATEMENTS 

REAL  DATA(0;  1000,7) 
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INTEGER  N 

*  PRINT  BANNER  AND  START  SUBROUTINE 

PRINT* 

PRINT*,1*****  RUNNING  ANALOG  DATA  FILE  SUBROUTINE  *****' 
PRINT* 

DO  100  I  =  1,  N 
DATA(I,6)  =  (N-I+1)*DATA(I,4) 
DATA(I,7)  =  (N-I)*DATA(I,4) 
100  CONTINUE 
END 

*  SUBROUTINE  DISDAT  * 

*  * 

*  SUBROUTINE  DISDAT  CREATES  LOAD  AND  DISPLACEMENT  BASED  ON  THE 

*  DATA  BUT  THE  DISPLACEMENT  IS  CONTROLLED  BY  THE  MACHINE  CROSSHEA    * 

*  SPEED  AND  THE  DATA  SAMPLING  RATE.   A  RANDOM  NUMBER  GENERATOR  IS 

*  USED  TO  SOLVE  FOR  THE  LOAD  IF  THE  DISPLACEMENT  OCCURS  AT  A  BREAK 

*  POINT. 

*  -v 

SUBROUTINE  DISDAT( DDATA , DATA , N , E , L , CHS , SR , POINTS ) 

*  DECLARATION  STATEMENTS 

REAL  DDATA(0;5000,2),DATA(0; 1000, 7) ,E, CHS ,SR,DELD, DISP, LIMIT, 
1  RANGE, RNG,L 
INTEGER  SEED  1, BREAKS, POINTS, N 

*  PRINT  BANNER  AND  START  SUBROUTINE 

PRINT* 

PRINT*,******  RUNNING  DISCRETE  DATA  SUBROUTINE  *****' 

PRINT* 

PRINT*, 'INPUT  ANY  ODD  INTEGER  (RNG  SEED)' 

READ*, SEED 1 

DELD  =  CHS/(SR*60) 

BREAKS  =  0 

*  SET  LIMIT  EQUAL  TO  DISPLACEMET  OF  FIRST  BREAK 

LIMIT  =  DATA(1,5) 
DO  100  I  =  1,5000 
DISP  =  DELD* I 
150  CONTINUE 

IF(DISP  . LT.  LIMIT)  THEN 

DDATA(I,1)  =  DISP 

DDATA(I,2)  =  ((N  -  BREAKS)*E*DISP)/L 

GOTO  199 

ENDIF 

IF(DISP  . EQ.  LIMIT)  THEN 

RANGE  =  (((N  -  BREAKS)*E*DISP)-((N-BREAKS-1)*E*DISP))/L 

DDATA(I,2)  =  (RNG(SEED1)*RANGE)+(((N-BREAKS-1)*E*DISP)/L) 

DDATA(I,1)  =  DISP 

GOTO  199 

ENDIF 

BREAKS  =  BREAKS  +  1 

IF( BREAKS  .EQ.  N)  GOTO  200 
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*  INCREMENT  TO  NEXT  BREAK  DISPLACEMENT 

LIMIT  =  DATA(BREAKS+1,5) 
GOTO  150 

199  CONTINUE 
100  CONTINUE 

200  CONTINUE 

*  POINTS  =  NUMBER  OF  DISCRETE  DATA  POINTS 

POINTS  =  I 
END 

*  SUBROUTINE  NOISE  * 

*  SUBROUTINE  NOISE  ALLOWS  THE  USER  TO  SIMULATE  MACHINE  NOISE  IN  THE  * 

*  DISCRETE  LOAD  DATA.   THE  USER  INPUTS  THE  RANGE  OF  THE  NOISE  AND  * 

*  A  RANDOM  NUMBER  GENERATOR  IS  USED  TO  FIND  A  NOISE  FOR  EACH  DATA  * 

*  LOAD  IN  BETWEEN  THE  GIVEN  RANGE. 

it  Vc 

SUBROUTINE  NOISE(DDATA,  POINTS) 

*  DECLARATION  STATEMENTS 

REAL  DDATA(0;5000,2),PH,PL,RNG 
INTEGER  I, Q,SEED2, POINTS 

*  PRINT  BANNER  AND  START  SUBROUTINE 

PRINT* 

PRINT*, ' *****   RUNNING  NOISE  SUBROUTINE  ****** 

PRINT* 

PRINT*, 'DO  YOU  WANT  TO  ENTER  NOISE  INTO  DISCRETE  DATA' 

PRINT*, 'IF  YES  TYPE  1,  IF  NO  TYPE  0' 

READ* , Q 

IF(Q  . NE.  l)GOTO  200 

PRINT*, 'INPUT  ANY  ODD  INTEGER  (RNG  SEED)' 

PRINT* 

READ*,SEED2 

PRINT*, 'INPUT  HIGH  NOISE  RANGE  AND  LOW  NOISE  RANGE   (GMS)' 

READ*,  PH,  PL 

DO  100  I  =  0,  POINTS 

DDATA(I,2)  =  (DDATA(I,2)  -( (PH-PL) /2) )  +  RNG( SEED2)*(PH-PL) 
100  CONTINUE 
200  CONTINUE 

END 

*  SUBROUTINE  OUTPUT  * 

*  SUBROUTINE  OUTPUT  WRITES  THE  ANALOG  DATA  AND  THE  DISCRETE  DATA     * 

*  ONTO  DRIVE -B  OUTPUT  FILES  * 

SUBROUTINE  OUTPUT( DATA , DDATA , N , E , L , ALPHA , BETA , POINTS , ADATA ) 
REAL  DATA(0; 1000 , 7) ,DDATA(0;  5000,2) ,E, ALPHA, BETA, L, 
1ADATA(0; 2000,2) 
INTEGER  N,P0INTS,Q1,Q2,I,II 

PRINT*,'*****  RUNNING  OUTPUT  SUBROUTINE  *****' 
PRINT* 
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PRINT*,' DO  YOU  DESIRE  A  LISTING  OF  THE  ANALOG  DATA' 

PRINT*,' IF  YES  TYPE  1,  IF  NO  TYPE  ZERO' 

READ*,Q1 

PRINT*, 'DO  YOU  DESIRE  A  LISTING  OF  THE  DISCRETE  DATA' 

PRINT*, 'IF  YES  TYPE  1,  IF  NO  TYPE  ZERO' 

READ* , Q2 

IF(Q1  .EQ.  1)  THEN 

*  WRITE  THE  ANALOG  DATA 

PRINT* 

OPEN(15,  FILE='ADATA') 

WRITE(15,  960) 

960  FORMAT(   /   '    DISP    LOAD') 
DO  105  I  =  1,N 

II  =  (2*I)-1 
ADATA(II,1)  =  DATA(I,5) 
ADATA(II+1,1)  =  DATA(I,5) 
ADATA(II,2)  =  DATA(I,6) 
ADATA( 11+1,2)  =  DATA(I,7) 

105  CONTINUE 

DO  106  I  =  0,2*N 

WRITE(15,961)  ADATA(I,1),  ADATA(I,2) 

106  CONTINUE 

961  FORMATC   F10. 6,2X,F10. 4) 
CLOSE(15) 

PRINT* 

PRINT*,'  ANALOG  DATA  IS  IN  FILE  "ADATA  FILE'" 

PRINT* 

END  IF 

IF(Q2  .EQ.  1)  THEN 

*  WRITE  THE  DISCRETE  DATA 

OPEN(15,  FILE='DDATA') 

WRITE( 15,951) 
951  FORMAT(/  'POINT       DISP        LOAD') 

DO  110  I  =  0, POINTS 

WRITE( 15 ,901) I ,  DDATA( 1,1) ,DDATA( 1,2) 
110  CONTINUE 

PRINT*,'  DISCRETE  DATA  IS  IN  FILE  "DDATA  FILE'" 

PRINT* 

CLOSE(15) 

END  IF 
901  FORMAT(   I5,2X,F10.6,  4X,F10.4) 

END 

FUNCTION  RNGCSEED)   RANDOM  NUMBER  GENERATOR  * 

FUNCTION  RNG(K) 

INTEGER  K,M,C0NST1 

REAL  RNG,  CONST2 

PARAMETERCC0NST1=2147483647,  CONST2=. 4656613E-9) 

SAVE 

DATA  M  /0/ 

IF(M  .EQ.  0)  M  =  K 

M  =  M*65539 

IFCM  .LT.  0)  M  =  (M+l)  +  CONST1 
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RNG  =  M  *   C0NST2 

END 
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APPENDIX  C 
COMPLIANCE    REMOVAL 

The  objective  of  the  analysis  of  fiber  bundle  testing  data  is  to 
determine  the  fiber  strength  distribution  parameters.  To  this  end,  the 
actual  fiber  displacement  data  must  be  obtained  from  the  total  test 
displacement  data.  A  significant  portion  of  the  test  displacement  data 
is  due  to  machine  compliance.  This  appendix  describes  a  method  to 
determine  this  displacement  and  subtract  it  from  the  test  displace- 
ment data  consisting  of  the  load  train  displacement  and  the  sample 
displacement. 

A.     ISOLATION   OF  COMPLIANCE 

Fiber  bundle  testing  is  accomplished  with  the  use  of  the  Instron 
Universal  Testing  Instrument  Model  4206  (INSTRON)  machine.  The 
fiber  test  records  displacement  and  the  load  measured  by  the 
INSTRON  load  cell.  The  test  is  controlled  by  software  that  operates 
the  machine  on  a  controlled  cross-head  speed  and  takes  data  at  a 
user-selected  sampling  rate.  The  machine  is  schematically  illustrated 
in  Figure  43. 

A  typical  bundle  test  load  vs  displacement  graph  is  shown  in  Fig- 
ure 44.  The  curve  may  be  partitioned  into  three  regions.  Region  I  is 
the  concave  non-linear  region.  This  is  concavity  is  due  to  the  compli- 
ance on  the  load  train  and  the  variation  in  gauge  length  (slack)  amount 
of  each  filament. 
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Figure  43 
INSTRON  Testing  Machine 
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Figure  44 
Typical  Bundle  Test  Load  vs.  Displacement  Graph 
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This  is  due  to  each  filament  of  the  bundle  not  having  exactly  the 
same  gauge  length;  the  curve  becomes  progressively  steeper  as  each 
filament  starts  to  carry  the  load.  Region  II  is  the  linear  region.  Here  all 
of  the  filaments  are  carrying  the  load  equally.  The  bundle  displace- 
ment is  a  linear  function  of  displacement  due  to  Hooke's  Linear 
Stress-Strain  law  shown  in  Equation  23. 

P  =  nE5 

5  =  P/nE  (23) 

where     P  =  load 

n  =  number  of  filaments 
E  =  filament  modulus 
8  =  displacement 

Region  III  is  the  failure  region  where  the  load  fluctuates  as  each 
fiber  breaks. 

The  load  data  that  the  INSTRON  records  is  derived  from  an  elec- 
trical signal  from  the  load  cell.  This  load  is  the  load  taken  by  the 
fibers,  but  also  includes  the  load  taken  up  by  the  testing  machine  load 
train.  This  load  causes  deformations  in  the  load  train  consisting  of  the 
gears  of  the  machine,  the  grips,  grip  base,  and  all  other  structural 
components.  These  components  of  the  load  train  combined  to  form 
machine  compliance.  The  INSTRON  shown  in  Figure  43  may  now  be 
modeled  as  shown  in  Figure  45.  The  INSTRON  is  equivalent  to  two 
springs  in  series,  the  machine  compliance  and  the  fibers. 
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Figure  45 
INSTRON  Testing  Machine  Equivalent  Model 

With  the  spring  in  series  model,  the  total  displacement  of  the 
system  due  to  a  load  is  equal  to  Equation  24.  The  load  due  to  the  fibers 
is  Equation  25  and  the  load  due  to  compliance  is  equation  26. 


8t  =  5f  +  5C  =  P/(Kf  +  Kc) 


(24) 


5f  =  P/Kf 


(25) 


6c  =  P/Kc 


(26) 


where     5t  =  total  bundle  displacement 

8f  =  fiber  displacement 

8C  =  compliance  displacement 

Kf  =  fiber  stiffness  coefficient 

Kc  =  compliance  stiffness  coefficient 
Equations  25  and  26  demonstrate  that  the  fiber  and  compliance 
displacements  are  functions  of  load.  If  the  stiffness  coefficients  and 
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the  load  are  known,  then  the  displacement  due  to  the  fibers  may  be 
found  by  Equation  27.  It  is  important  to  note  that  the  stiffness  coeffi- 
cients for  the  fibers  and  the  compliance  are  likely  to  be  nonlinear.  The 
fiber  nonlinearity  is  primarily  due  to  the  slack  in  the  early  displace- 
ment. The  compliance  will  be  nonlinear  due  to  the  physical  processes 
that  contribute  to  form  the  compliance.  Many  of  these  processes  are 
directional  in  nature. 


8(P)f  =  8(P)t  -  5(P)C 


(27) 


From  Equation  27,  the  typical  bundle  test  graph  in  Figure  1  may 
now  be  represented  as  Figure  46  with  the  compliance  displacement 
being  subtracted  from  the  fiber  displacement. 
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Figure  46 

Typical  Bundle  Test  Load  vs  Displacement  Graph  Compliance 
and  Fiber  Displacement  Separated 
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To  determine  the  fiber  displacement,  a  means  to  subtract  the 
compliance  displacement  is  required.  An  equation  must  be  deter- 
mined to  fit  8(P)C  to  implement  Equation  27. 

If  Kc  were  linear,  it  would  be  a  simple  matter  to  solve  for  the  8(P)C 
function.  However,  as  mentioned  above,  the  function  will  not  be  linear. 
Because  of  this,  8(P)C  must  be  measured  by  experimentation.  The  data 
must  then  be  curve-fit  to  form  a  continuous  function;  this  function  may 
then  be  used  to  implement  Equation  27. 

B.     COMPLIANCE   TESTING 

From  the  model  illustrated  in  Figure  45,  a  means  to  eliminate  the 
displacement  of  the  fibers  would  be  to  set  the  gauge  length  of  the 
bundle  to  zero.  The  fibers  may  not  be  eliminated  because  the  gripping 
mechanism  to  hold  the  bundle  in  place  is  a  part  of  the  compliance,  but 
the  gauge  length  can  be  brought  very  nearly  to  zero.  To  do  this,  a  sam- 
ple is  prepared  like  the  one  shown  in  Figure  47.  If  there  are  no  fibers 
that  displace,  all  of  the  displacement  recorded  must  be  due  to 
machine  compliance. 

To  measure  5(P)C,  a  bundle  test  is  completed  with  the  0  gauge 
length  sample  in  accordance  with  the  procedures  outlined  in 
Appendix  D.  A  typical  INSTRON  compliance  test  load  vs.  displacement 
graph  is  shown  as  Figure  48. 
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Figure  48 
INSTRON  Compliance  Test  Load  vs.  Displacement  Graph 
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The  first  region,  0  to  Del,  is  very  nonlinear.  This  is  due  to  the 
directional  compliance  in  the  INSTRON  as  the  displacement  is  first 
applied.  Physical  processes  include  the  backlash  in  the  gears,  nonlin- 
ear friction  and  damping  as  the  gears  and  clutches  engage,  and  noise 
below  the  sensitivity  of  the  load  cell. 

Region  two,  Del  to  Dc2,  is  the  linear  region.  The  compliance  in 
this  region  is  due  to  linear  spring  tension  of  the  machine  as  modeled 
in  Figure  44.  The  K  is  very  high,  which  is  expected  for  a  testing 
machine.  It  should  be  high  so  that  the  load  it  absorbs  is  very  low  when 
compared  to  the  load  being  absorbed  by  the  sample. 

Region  three  is  nonlinear.  It  is  believed  that  in  this  region  the 
grips  are  beginning  to  fail  and  that  the  fibers  are  actually  slipping 
between  the  copper  tabs.  The  data  in  this  region  does  not  reflect  the 
machine  compliance  and  is  not  utilized  in  the  analysis.  It  is  assumed, 
however,  that  if  the  sample  had  not  slipped  the  compliance  curve 
would  have  continued,  extending  the  linear  region. 

C.     CURVE-FITTING  PROCEDURE 

The  fitting  of  a  curve  to  approximate  P(8)c  was  completed  by 
separating  the  function  into  two  regions,  selecting  a  power  series 
model,  Equation  28,  for  region  I,  a  linear  model;  Equation  29,  for  sec- 
tion II;  and  setting  boundary  conditions  at  Del,  Equations  30  and  31. 

forO<Dc<Dcl  Dc  =  Al(P)n  +  A2  (28) 

forDcl<Dc  Dc=A3P  +  A4  (29) 
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Boundary  Conditions 

(1)  Del  (Equation  28)  =  Del  (Equation  29)  (30) 

(2)  Del'  (Equation  28)  =  Del'  (Equation  29)  (31) 
If  Dc2  and  Del  are  known,  then  from  Equation  29  for  region  two: 

A3  =  Dc2  -  Dcl/Pi  -  p(Dcl)  (32) 

A4  =  -A3P1  +  Del  (33) 

From  Equations  28,  29,  and  31,  an  equation  for  Ai  can  be  derived: 

3  Dc/3  P  =  nAi(P)n-i 

a  Dc/a  p  =  A3 

A3  =  nAi(P)n-! 
Ai  =A3/n  (Pi)n-l  (34) 

From  Equations  29,  30,  and  34,  an  equation  for  A2  can  be  derived: 

Del  =  ((A3/n  (Del  -  DcO)11'1)  *  (Del  -  Dc0)n/A3  +  A2 

A2  =  Dci -A3Pi/n  (35) 

This  gives  us  five  unknowns  and  four  equations.  Equations  32  and 
33  can  be  solved  directly  from  information  gained  from  the  0  gauge 
length  data.  Equations  34  and  35  are  coupled  by  n,  so  A3  and  A4  may 
be  solved  for  by  optimizing  the  variable  n.  To  optimize  n,  an  error 
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function  is  created  and  its  derivative  is  taken  and  then  summed  for  all 
the  data  points  between  zero  and  Del  for  different  values  of  n.  When 
the  error  function  is  minimized  to  approximately  zero,  n  is  known. 
The  n  can  then  be  substituted  into  equations  13  and  14  to  yield  the 
constants  A3  and  A4. 

Substituting  Equations  32  and  33  into  Equation  28  gives  Equation 
36. 

Dc(P)  =  (A3P1/11)  *  ((P/Pi)11-1)  +  Del  (36) 

3  Dc(Pi)/3  N  =  (A3P1/11)  *  ln(P/Pi)  *  (P/Pi)n 

+  AaPin-2(l  -  (P/Pjn)  (37) 

The  error  function  becomes 

H(n)  =  I  (Dc(Pi)  -  Di)  9  Dc(Pi)/a  N  (38) 

for  i  =  1,2,3  ...  number  of  data  points  from  0  to  Del.  Equation  38  is 
then  iterated  until  H(n)  goes  to  zero  (.0000). 

The  process  of  solving  for  the  coefficients  is  greatly  enhanced  by 
the  use  of  spreadsheets  and  microcomputers.  Programs  such  as  Lotus 
1-2-3  and  Microsoft  Excel  allow  the  user  to  import  the  data  directly 
into  a  spreadsheet  and  manipulate  the  data.  The  results  for  the 
INSTRON  compliance  used  in  this  project  were  found  using  Lotus 
1-2-3. 
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A  Fortran  program  is  in  development  to  solve  this  procedure  but 
is  not  functional  at  this  time. 
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APPENDIX  D 
BUNDLE  TESTING 

The  failure  load  for  graphite  bundles  was  determined  with  the  use 
of  the  Instron  Universal  Testing  Instrument  Model  4206  (INSTRON) 
and  associated  IBM  INSTRON  control  software,  version  4.01.  The  fiber 
strength  statistics  were  found  by  removing  the  compliance  displace- 
ment with  the  software  developed  in  this  appendix.  The  compliance 
displacement  function  was  approximated  by  the  curve  fit  detailed  in 
Appendix  C.  The  modified  data  was  then  used  with  software  developed 
by  Lt.  Joseph  Schmidt  to  determine  the  fiber  parameters  a  and  p  for 
the  given  gauge  length. 

All  testing  was  completed  in  the  Advanced  Composites  Laboratory, 
Naval  Postgraduate  School,  Monterey.  Laboratory  Technician  Mr.  Jim 
Nageotte  and  Professor  Edward  Wu  assisted  with  the  preparation  of 
samples  and  the  testing. 

A.      BUNDLE  TESTING  PROCEDURE 

The  testing  was  completed  in  accordance  with  the  following 
procedure. 

•  Turn  on  INSTRON  and  allow  60  minutes  to  warm  up. 

•  Calibrate  the  INSTRON  in  accordance  with  posted  procedures  for 
mechanical  and  electrical  calibration  (only  required  for  first  test 
of  the  day). 

•  Turn  on  attached  IBM  AT  and  enter  INSTRON  test  software.  From 
C  prompt,  type:  MT\DATA\MT  <enter> 

•  Load  sample  into  machine. 
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•  Ensure  sample  is  vertically  aligned  correctly;  adjust  if  necessary. 

•  Set  zero  load  level. 

Depress  Load  Bal  button  and  then  enter  on  the  INSTRON. 

•  Set  cross-head  speed  to  .2  mm/min. 

•  Set  gauge  length. 

Run  cross-head  up  until  a  load  of  .2  ±.02  kg  is  achieved. 
Depress  Gauge  Length  button  on  the  INSTRON. 
Run  cross-head  down  until  a  load  of  zero  is  achieved. 
Run  cross-head  up  until  a  load  of  .2  ±.02  kg  is  achieved. 
Extension  should  read  0.0.  If  it  does  not,  repeat  all  steps. 
Run  cross-head  down  until  a  load  of  zero  is  achieved. 
Depress  jog  down  switch  twice. 
Depress  Gauge  Length  button. 
*(Crosshead  Speed  may  be  adjusted  if  necessary) 

Select  appropriate  test  method. 

Enable  IEEE  interface. 

Depress  IEEE  button  on  the  INSTRON. 

Run  the  test. 

Plot  the  results. 

Plot  on  HP  plotter  using  Plot  option  in  INSTRON  software. 

Put  data  into  ASCII  file. 

Dump  data  to  hard  disk  or  floppy  disc  using  Utilities  option  in 
INSTRON  software. 

Exit  INSTRON  software. 

Remove  Compliance  displacement. 

Rename  test  data  file  to  TEST. DAT  and  put  onto  floppy  disk. 

Place  data  disk  in  drive  B: 
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•  Place  compliance  removal  program  disk  in  drive  A: 

•  Type  T5  <enter> 

The  disk  in  drive  B:  now  contains  the  data  with  compliance 
removed  and  the  data  with  compliance  not  removed.  The  data  is  now 
ready  to  be  run  in  the  software  developed  by  Lt.  Joseph  Schmidt  that 
will  retrieve  the  fiber  parameters  a  and  p. 

B.     PROGRAM  LISTING 

The  compliance  removal  program  "program  T5"  takes  the  five 
coefficients  for  the  curve  fit  compliance  function  and  then  for  the  cor- 
responding load  subtracts  the  displacement  from  the  data  due  to 
compliance.  The  method  for  determining  the  compliance  curve  fit 
coefficients  is  detailed  in  Appendix  C. 

The  program  is  written  in  Microsoft  Fortran  4.01.  It  will  run  on 
any  IBM/compatible  personal  computer.  To  operate  the  program, 
place  the  program  disk  into  the  A:  drive.  The  data  disk  must  be  placed 
into  the  B:  drive.  When  this  is  completed,  type:  T5  <enter>. 

PROGRAM  T5 

SUBROUTINE  LISTING: 

I  NIT  Initializes  arrays  to  zero 

LOAD  Opens  the  data  file  and  reads  in  the  machine  data 

CONVERT  Converts  machine  data  into  SI  units 

REMCON  •  Prompts  the  user  to  input  the  compliance  curve 

coefficients 

•  Removes  the  displacement  from  the  data  due  to 
compliance 
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OUTPUT 

VARIABLE  LISTING: 

DDISP(6000) 

DLOAD(6000 

INDEX 

XCONV 

YCONV 

Al 

A2 

N 

A3 

A4 

DC1 

ZERO 


Writes  to  the  B:  drive  the  modified  data 

Test  displacement  data 

Test  load  data 

Test  data  point  counter 

Displacement  conversion  factor 

Load  conversion  factor 

Compliance  curve  coefficient 

Compliance  curve  coefficient 

Compliance  curve  coefficient 

Compliance  curve  coefficient 

Compliance  curve  coefficient 

Transition   displacement  between   two   curve   fit 
equations 

Noise  level  of  load  cell 
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PROGRAM  T5 

VoWoWr  Vr  Vr  V- Vc  Vw  r  Vc  Vr -/w  c  Vr  Vr  Vc  VoV  VoV  V"V  Vc  ^ 

*  PROGRAM  T5  * 

*  * 

THIS  PROGRAM  WILL  TAKE  INSTRON  4200  SERIES  BUNDLE  RAW  DATA  OUTPUT  * 
AND  ALLOW  THE  USER  TO  CONVERT  THE  DATA  INTO  DESIRED  UNITS  AND  TO  * 
TO  REMOVE  MACHINE  COMPLIANCE  DISPLACEMENT.  THE  INPUT  DATA  FILE  * 
MUST  BE  PLACED  IN  THE  B-  DRIVE.   THE  INPUT  DATA  FILE  MUST  BE 

*  NAMED  TEST. DAT.   A  MAXIMUM  OF  6000  DATA  POINTS  IS  ALLOWED.   TO 

*  EXECUTE  THE  PROGRAM,  PLACE  THE  PROGRAM  DISK  IN  THE  A  DRIVE,  * 
TYPE  A-T5  AND  HIT  ENTER.  tHE  USER  SHOULD  THEN  FOLLOW  THE  SCREEN  * 
PROMPTS. 

THE  REQUIRED  INPUTS  ARE- 

*  -   NUMBER  OF  DATA  POINTS 

-  CONVERSION  FACTOR  FOR  LOAD  AND  DISPLACEMENT 

-  COMPLIANCE  CURVE  PARAMATERS  (IF  DESIRED)  * 

*  * 

THE  OUTPUT  DATA  FILE  WILL  BE  PUT  ON  THE  DISK  IN  THE  B  DRIVE  AS 

*  FILE  EXPER. OUT.    A  SUMMARY  PRINTOUT  OF  THE  INPUT  IS  PLACED  ON 

THE  DISK  IN  THE  B  DRIVE  AS  DATCON.  PRT  .  * 

*  * 

REAL  XCONV,YCONV,A1,A2,A3,A4,N,DDISP(6000),DLOAD(6000), 
+ZER0 
INTEGER  I,J,DP,INDEX(6000),Q 
PRINT*,1***  RUNNING  PROGRAM  T5  ***' 
PRINT* 

*  CALL  SUBROUTINE  INIT  TO  INITIALIZE  ARRAYS 

CALL  INIT(DDISP,DLOAD, INDEX, XC0NV,YC0NV) 

*  CALL  SUBROUTINE  LOAD  TO  LOAD  IN  DATA  INTO  PROGRAM 

CALL  LOAD(DDISP,DLOAD, INDEX, DP) 

*  ALLOW  USER  TO  CO! "VERT  DATA  INTO  DESIRED  UNITS 

PRINTS,'     DO  YOU  NEED  TO  CONVERT  DATA  INTO  [SI]  UNITS' 

PRINT*,'  THIS  SHOULD  ONLY  BE  REQUIRED  FOR  MACHINE  TEST  DATA' 

PRINT*,'         IF  YES  TYPE  1,  ELSE  TYPE  0  TO  CONTINUE' 

READ*,Q 

IF(Q  . EQ.  1)  THEN 

*  IF  USER  DESIRES  TO  CONVERT  UNITS  CALL  SUBROUTINE  CONVERT  TO  CONVERT 

*  DATA 

CALL  CONVERT(XCONV,YCONV,DP,DDISP,DLOAD) 
END  IF 

*  CALL  SUBROUTINE  REMCOM  TO  REMOVE  COMPLIANCE  FROM  DATA 

CALL  REMCOMC DDI SP , DLOAD , A 1 , A2 , A3 , A4 , N , DP , ZERO , DC  1 ) 

*  CALL  SUBROUTINE  OUTPUT  TO  OUTPUT  DATA  INTO  DATA  FILES  ON  DISK  B 

CALL  OUTPUTC INDEX, DD ISP, DLOAD, A 1,A2, A3, A4,N, DP, ZERO, XCONV,YCONV, 
+DC1) 
PRINT*,'***  END  OF  PROGRAM  ***' 
END 

*  SUBROUTINE  INIT 

A  * 

*  THIS  PROGRAM  WILL  INITIALIZE  THE  DATA  ARRAYS,  INDEX  AND  CONVERSION  * 

*  FACTORS  PRIOR  TO  INPUTING  NEW  DATA.  * 

SUBROUTINE  INIT(DDISP, DLOAD, INDEX, XCONV,YCONV) 
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REAL  DDISP( 6000), DLOAD( 6000), XCONV, YCONV 
INTEGER  INDEX(6000) 

*  PRINT  HEADER 

PRINT''",  '***  RUNNING  SUBROUTINE  INIT  ***' 
PRINT* 

*  RUN  LOOP  TO  SET  DATA  ARRAY  INDEX,  DDISP  AND  DLOAD  TO  0 

DO  10  I  =  1,6000 
INDEX(I)  =  0 
DDISP(I)  =  0. 
DLOAD(I)  =  0. 
10  CONTINUE 

*  SET  CONVERSION  FACTORS  TO  1.0 

XCONV  =  1. 
YCONV  =  1. 
END 

*  SUBROUTINE  LOAD 

*  -V 

SUBROUTINE  LOAD  WILL  OPEN  DATA  FILE  B -TEST. DAT  AND  READ  IN  THE  * 

*  DATA  WITH  AN  UNFORMATTED  READ  STATEMENT.  THE  NUMBER  OF  DATA  * 
POINTS  (DP)  TO  BE  READ  MUST  BE  ENTERED  CORRECTLY  OR  PROGRAM  WILL  * 
ERROR  OUT  IN  THIS  SUBROUTINE. 

&;fc*&sW«V}VsV:toWfAAA»V>V»V}V3YA*iWc»V&ywrV«V 

SUBROUTINE  LOAD(DDISP, DLOAD , INDEX, DP) 
REAL  DDISP(6000),DLOAD(6000) 
INTEGER  I,DP,INDEX(6000) 

*  PRINT  HEADER 

PRINT*,'***  RUNNING  SUBROUTINE  LOAD  ***' 
PRINT* 

*  OPEN  DEVICE  15,  FILE  B-TEST. DAT 

OPEN( 15 ,FILE=' B-TEST. DAT' ) 

*  INPUT  NUMBER  OF  DATA  POINTS  (DP) 

PRINT*, 'INPUT  NUMBER  OF  DATA  ELEMENTS  TO  BE  READ' 
PRINT*,'   @@@  MAXIMUM  INPUT  VALUE  IS  6000  @@@' 
RE AD*, DP 

*  RUN  LOOP  TO  READ  IN  DATA  POINTS  INDEX,  DDISP  AND  LOAD 

DO  10  I  =  1,DP 

READ(  15  ,  ••')  INDEXC  I)  ,DDISP(  I)  ,DL0AD(  I) 
10  CONTINUE 

*  CLOSE  DEVICE  15 

CLOSE  (15) 
END 

*  SUBROUTINE  CONVERT  * 

*  * 

*  SUBROUTINE  CONVERT  WILL  CONVERT  THE  DATA  FROM  FILE  TEST.  DAT  INTO 

*  WHICHEVER  UNITS  THE  USER  DESIRES.   THE  USER  MUST  INPUT  CONVERSION  * 

*  FACTORS  WHICH  WILL  BE  DIVIDED  FROM  THE  ORIGINAL  DATA.   THE  X  CONV  * 

*  CORRESPONDS  TO  DISPLACEMENT  AND  THE  Y  CONV  CORRESPONDS  TO  LOAD  OR  * 

*  FORCE.   IF  NO  CONVERSION  IS  DESIRED,  ENTER  1  FOR  THE  CONVERSION 

*  FACTORS.  * 

SUBROUTINE  CONVERT( XCONV , YCONV , DP , DDI SP , DLOAD ) 
REAL  DDISP( 6000), DLOAD( 6000), XCONV, YCONV 
INTEGER  I, DP 

*  PRINT  HEADER 
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PRINT*,'***  RUNNING  SUBROUTINE  CONV  ***' 
PRINT* 

PRINT--'",  'PROGRAM  T5  WILL  CONVERT  DATA  INTO  [SI]  UNITS' 
PRINT* 

*  INPUT  X  OR  DISPLACEMENT  CONVERSION  FACTOR 

PRINT*, 'INPUT  X  CONVERSION   (DISPLACEMENT)" 
READ*,  XCONV 

*  INPUT  Y  OR  LOAD/FORCE  CONVERSION  FACTOR 

PRINT*, 'INPUT  Y  CONVERSION   (LOAD/FORCE)' 
READ*,  YCONV 

*  CONVERT  DATA 

DO  10  I  =  1,DP 
DDISP(I)  =  DDISP(I)/XCONV 
DLOAD(I)  =  DLOAD(I)/YCONV 
10  CONTINUE 
END 

*  SUBROUTINE  REMCON 

-v  * 

*  SUBROUTINE  REMCON  WILL  IF  THE  USER  DESIRES  REMOVE  THE  MACHINE 

*  COMPLIANCE  DISPLACEMENT  FROM  THE  DATA.   THE  USER  MUST  INPUT  THE 

*  COMPLIANCE  CURVE  PARAMETERS  AND  THE  DATA  ZERO  LOAD  LEVEL.  * 

*  THE  USER  SHOULD  TAKE  CARE  TO  ENSURE  THE  COMPLIANCE  PARAMETERS  ARE  * 

*  THE  SAME  UNITS  AS  THE  CONVERSION  FACTORS. 

SUBROUTINE  REMCOM(DDISP,DLOAD,Al ,A2 , A3 ,A4,N,DP,ZER0,DC1) 
REAL  DDISP(6000),DLOAD(6000),A1,A2,A3,A4,N,ZERO,DC1 
INTEGER  I,DP,Q 

*  PRINT  HEADER 

PRINT*,'***  RUNNING  SUBROUTINE  REMCOM  ***' 
PRINT* 

*  ALLOW  USER  TO  REMOVE  COMPLIANCE  FROM  THE  DATA 

PRINT*, 'DO  YOU  DESIRE  TO  REMOVE  THE  MACHINE  COMPLIANCE  FROM  THE  DA 
+TA' 
PRINT*, 'IF  YES  TYPE  1,  ELSE  TYPE  0  TO  CONTINUE' 
READ*,Q 
IF(Q  . EQ.  1)  THEN 

*  IF  USER  DESIRES  TO  REMOVE  COMPLIANCE,  INPUT  COMPLIANCE  PARAMETERS 

PRINT*, 'INPUT  COMPLIANCE  CURVE  COEFFICIENTS' 

PRINT*, 'INPUT  Al' 

RE AD*, A 1 

PRINT*,' INPUT  A2' 

READ*,A2 

PRINT*, 'INPUT  N' 

READ* , N 

PRINT*,* INPUT  A3' 

RE AD*, A 3 

PRINT*, 'INPUT  A4' 

READ* , A4 

PRINT*, 'INPUT  DC1' 

RE AD*, DC  1 

PRINT*,' INPUT  ZERO  NOISE  LEVEL  (LOAD  IN  KG)' 

READ*, ZERO 

*  REMOVE  COMPLIANCE  DISPLACEMENT  FROM  THE  DATA 

PRINT* 

PRINT*,'***  REMOVING  MACHINE  COMPLIANCE  FROM  DATA  ***' 
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PRINT* 

DO  100  I  =  1,DP 

*  DETERMINE  IF  DISPLACEMENT  DATA  IS  IN  COMPLIANCE  POWER  CURVE  OR 

*  LINEAR  REGION  AND  SUBRTRACT  CORRECT  AMOUNT  OF  DISPLACEMENT 

IF(DLOAD(I)  .  LT.  ZERO)  THEN 
GOTO  99 
END  IF 

IF(DDISP(I)  .  LT.  DC1)  THEN 

DDISP(I)  =  DDISP(I)  -  ((Al*(DLOAD(I)**N))  +  A2) 
ELSE 

DDISP(I)  =  DDISP(I)  -  ((A3*DLOAD(I))  +  A4) 
ENDIF 
99  CONTINUE 
100  CONTINUE 
ENDIF 
CONTINUE 
END 

*  SUBROUTINE  OUTPUT  * 

*  SUBROUTINE  OUTPUT  WILL  WRITE  THE  DATA  TO  FILE  EXPER. OUT.   IT  WILL  * 

*  ALSO  CREATE  A  DATA  SUMMARY  FILE  AS  DATCON. PRT  .   THE  OUTPUT  FILES  * 

*  ARE  PLACED  ONTO  THE  DISK  IN  THE  B  DRIVE. 

SUBROUTINE  OUTPUT( INDEX , DDI SP , DLOAD , Al , A2 , A3 , A4 , N , DP , ZERO , XCONV , 
+YC0NV,DC1) 
REAL  DDISP(6000),DLOAD(6000),A1,A2, A3, A4,N, ZERO, XCONV, YCONV 
INTEGER  I,DP,INDEX(6000),Q 

*  PRINT  HEADER 

PRINT*,1***  RUNNING  SUBROUTINE  OUTPUT  ***' 
PRINT* 
10   PRINT*,' DATA  WILL  BE  OUTPUT  TO  FILE  B-EXPER.  OUT' 

*  PROMPT  USER  TO  PLACE  TARGET  DATA  DISK  IN  DRIVE  B- 

PRINT*, 'ENSURE  A  DISK  WITH  SUFFICIENT  SPACE  IS  IN  DRIVE  B-' 

PRINT*, 'TYPE  1  TO  CONTINUE' 

READ* , Q 

IF(Q  .  NE. 1)  GOTO  10 

*  OPEN  DEVICE  12  AND  WRITE  DATA  TO  FILE  B-EXPER. OUT 

OPEN(12,  FILE=' B-EXPER.  OUT') 
DO  100  I  =  1,DP 

WRITE(12,900)INDEX(I),DDISP(I),DLOAD(I) 
100  CONTINUE 

*  CLOSE  DEVICE  12 

CLOSE(12) 

*  ALLOW  USER  TO  CREATE  A  DATA  SUMMARY  FILE 

PRINT*, 'DO  YOU  DESIRE  A  SUMMARY  PRINTOUT  FOR  THE  DATA  CONVERSION' 
PRINT*, 'IF  YES  TYPE  1,  ELSE  TYPE  0  TO  EXIT' 
READ* , Q 

*  IF  USER  DESIRES  DATA  SUMMARY  FILE,  OPEN  DEVICE  12  AND  WRITE 

*  PARAMETERS  TO  FILE  DATCON. PRT 

IF(Q  . EQ. 1)  THEN 

OPEN(12,  FILE='B -DATCON.  PRT') 

WRITEC 12,910) 

WRITE( 12, 911)DP, XCONV, YCONV 

WRITEC 12, 9 12) ZERO, DC 1 

WRITEC 12, 913)A1,A2,N, A3, A4 
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CL0SE(12) 
END  IF 
CONTINUE 
WRITE  FORMAT  STATEMENTS 
900  FORMAT(   16 , IX, F9. 4 , IX, F9.  4) 

"COMPLIANCE  REMOVAL  DATA  CONVERSION  COEFFICIENTS',//) 


910  FORMAT ( 

911  FORMAT( 
1 

2 

912  FORMAT ( 
1 

913  FORMAT( 
1 

2 
3 
4 
5 
6 
7 
END 


,8X,I7/ 
,5X,F10.  7/ 
,5X,F10. 7//) 
,10X,F7.4/ 
,8X,F9.4//) 


NUMBER  OF  DATA  POINTS 
X  CONVERSION  FACTOR 
Y  CONVERSION  FACTOR 
ZERO  NOISE  LEVEL  (KG) 
DC1  -LINEAR  SECTION- 
COMPLIANCE  CURVE  COEFFICIENTS'// 
POWER  FACTOR  PORTION'// 
Al  =    '  ,F6.4/ 
A2  =    '  ,F6.4/ 
N  =    '  ,F6.4/ 
LINEAR  PORTION'// 
A3  =    '  ,F6.4/ 
A4  =    '  ,F6.4/) 
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APPENDIX  E 
TEST  DATA 

A.      COMPLIANCE  TEST  DATA 

The  INSTRON  compliance  test  data  forms  are  contained  in  this 
appendix.  They  are  tabulated  in  Test  Number  sequence.  The  raw  data 
is  stored  on  the  IBM  AT  personal  computer  in  the  Composites  Labora- 
tory, Naval  Postgraduate  School,  in  director  MT/DATA/  as  file  (Test 
Number).  MAD. 
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NFS  Cooposi t es  Lab 

U.S. Naval  Postgraduate  School 
Monterey,  Cfl 

Fiber  Bundle  Compliance  Curve 

Zero  guage  length  ' 

. 15  as/tin 

Test  type:     Compliance  Curve  Instron  Corporation 

Series  11  Autoiated  Materials  Testing  System  v't.Ok 
Operator  name:  Jim  Test  Date:  August  5,1988 

Sample  Identification:  0GL1O50B  Sample  Type:   bndl 

Interface  Type:  4200  Series 

Machine  Farameters  of  test: 

Sample  Rate  (pts/sec):  10.00  Humidity  (  X  ):    40 

Crosshead  Speed  (in/min  ):    .150  Temperature  (deg.  F):    70 

Dimensions: 

Spec.  1 

Lin.  Density  (den)        8036.0 
Gauge  length  (ei)  1.0000 

Gut  of  1  specnens,  0  excluded. 

Load  Displcrent  X  Strain  Displcsent  Load 

Modulus      at  at        at        at  at 

Specimen    label1              Maximut  Break      Break     MaximuJ  Maximum 

Nu*ber               (gn/den)     (Kg)  (as)       (X)        (r.j)  (Kg) 

1      good      99.95      29.57      .3011      30.11      .89<i8      29.57 

Mean:  99.95  29.57  .3011  30.11  .2948  89.57 

Standard 

Deviation:  

Mean  - 

2.00  »  Sdv:  

Mean  + 

2.00  »  Sdv:  -- --     -- 

Mimrum:  99.95  29.57  .3011  30.11  .29«  29.57 

Maximum:  97.95      29.57      .3011      30.11      .8942      89.57 
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VERSION   't.Olr 
Samp  1  e  j  d 
Vcrr.  i  on 
Repor  t  f  i  1  ett 


OCL  1  Of;  00 
t\  .01c 
34 


Test  cJ  a  t  e 

Vers  i  on  date 

Operator 


August   5, 19QG 
06/09/BB 

Jim 


X  conver  r;  i  on 
Y  conversion 


.03937008 
2.20462300 


X  A/D  offset 
Y  A/D  offset 


.  O'.x.'i-i 
,  0<  *0<  > 


Samp  1  e  rate 

A/D  range 

Ca  1  lb  type 

Cal i ta  1 oad 

Tempera ture 

Test  type 

Bar  type 

Break  check 

Load  1  i  in  i  t 


1  0 .  00 
0 
ALII  DM  AT  I C 

1 10. 231 1 
70 
DENIER 
bndl 

1 00 . 00000 
1 10.231 10 


Ex  t ens o meter 

An  tos tar  t 

Geome  br  y 

Ca  ]  l  b  <"•::  tei  IS 

Hi  imi  ri  )  ty 

tt  spec i  mens 

Entry  cl  l  mens 

Thresh  de 1  ay 

1.::  tens  limit 


NO 

orr 

CYLINDRICAI 

.  000'  i 
60 

1 
YES 

.02205 
.  157'iH 


Samp  1 e  d  i mens i ons 
A :        . 0394    B : 


2036 . 0000 


,0394 


D: 


.039'+    E:  NO 


Spec i  men  0 

Max  j  mum  ]  oad 

Ma:;  extens 


65. 19  1 

.  (.)  1  3 


Stop  coll.  statu; 
Ha::  ]  oad  point  It 
Max  extens  pnt  tt 


DREA1  DETECT 
1  206 

1  309 


Number  of  elem e n t  s 


1309 


Specimen  dimensions  : 

A  :        .03 9 3 7 0 1  B :   2 0 3 6 . 0 0 0 00 00  C 


,0393701  D: 


,0393701  E:NL) 


Auxiliary  Input  Array  tt  1 
*******  *******  *******  *******  *******  *******  *******  *******  ******* 

Auxiliary  Input  Array  ti  2 
**.  *.***..$  *******  ********  *-«****-?-  *******  *******  *******  *******       .0  1 
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B.     BUNDLE  TEST  DATA 

The  INSTRON  bundle  test  data  forms  are  contained  in  this 
appendix.  They  are  tabulated  in  Test  Number  sequence.  The  raw  data 
is  stored  on  the  IBM  AT  personal  computer  in  the  Composites  Labora- 
tory, Naval  Postgraduate  School,  in  directory  MT\DATA\  as  file  (Test 
Number).MAD. 
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COHOSHES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939<.0 

GRAPHITE  BUNDLE  TEST  METHOD  170  CHI 

Test  lethod  for  Instron  ^206.  Used  for  testing  AS-<i 

Graphite  bundles  (appx  1000  fibers) 


Test  type: 


OIL  BUNDLE 


Operator  name:  MARK  JONES 


Instron  Corporation 

Series  II  Automated  Materials  Testing  Systen  v<i.0lc 

Test  Date:  September  B,  198B 


Sample  Identification:  080901 

Interface  Type:  <i200  Series 

Machine  Parameters  of  test: 

Satple  Rate  (pts/sec):   6.67 
Crosshead  Speed  lufim  ):   7.000 


Satple  Type:   AS  - ^ 


Humdity  (  X  ):        50 
Tetperature  (deg.  F):        73 


'BATCH  ID' 


0B0901  70  CM  WITH  OIL 


Spec.  1 


Lin.  Density  (tex) 
Gauge  length  (in) 


2036. 0 
700.00 


Out  of  1  specUens,  0  excluded. 


Bisplctent  X  Strain  Load 

at  at  at 

Speciien    Maxiaui  Maviiur  Maxitun 

Number       (n)  (X)  UN) 


1 


5.2^0 


.7<iB& 


,1359 


Mean: 


5.B«i0 


-7*186 


.1359 


Standard 
Deviation: 


117 


GRAPHITE  BUNDLE  TEST  (5  CM) 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939<t0 


Test  type:     Yarn/Fiber 
Operator  name:  Mark  Jones 

Sample  Identification:  090901 

Interface  Type:  'ic'OO  Series 

Machine  Paraieters  of  test: 

Sample  Rate  tpts/sec):   3.33 
Crosshead  Speed  (ni/iin  ):    .500 


Instron  Corporation 

Series  II  Automated  Materials  Testing  System  v't.Olc 

Test  Date:  September  9,1938 

Sample  Type:   AS-<i 


Humidity  (  J  ):    50 
Temperature  (deg.  P): 


73 


'BATCH  ID1         090901 
'OIL  OR  NO  OIL'     NO  OIL 

Dimensions: 

Spec,  I 

Lin.   Density   (tex)  2036.0 

Gauge  length   [nisi  50.000 

Out  of  1  specimens,  0  excluded. 

Displcment  X  Strain     Load 

at  at       at 

Specimen    Maximum  Maximum    Maximum 

Number       (m»)  (X)  .   IKN) 

1      .6500      1.300      ,23<t9 

Mean:      .6500      1.300      .E3<t9 

Standard 

Deviation:  

Mean  - 

2.00  i   Sdv:  -- - - 

Mean  ♦ 

8.00  »  Sdv:  

Minimus:    .6500      1.300      .83*19 

Maximus:         .6500  1.300  .53<i9 
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VERSION   4.01c 

Sample  id  :  090901 
Vers  i  an  :  4 .  0 1  c 

Report  filett  :  IE 


Test  date 

Vers i on  date 

Oper  a tor 


September   9, 1700 
06/ 09 /OH 
Marl:  Jones 


X  conversion  :     .03937008 
Y  conversion  :   E.E0462300 


X  A/D  offset 
Y  A /I)  offset 


,  IJl.K  M  1 

,  00i  K") 


Samp le  rate 

A /I)  range 

Ca 1 i  b  type 

C  a  1  i  b  1  o  a  d 

Tempers  lure 

Test  type 

B  a  r  t  y  p  e 

Break  check 

L.  o  a  d  limit 


3.33 
0 
AUTOMATIC 

1 10.E31 1 
73 
DENIER 

1 0 . 00000 
1  IE. 404 50 


E::  t ens o meter 

Autos  tart 

Geome t  ry 

Ca  1  j  b  e::  tens 

Humid  i  ty 

tt  spec  i  mens 

En t ry  d i mens 

Thresh  de ) ay 

E  >:  t  e  n  s  limit 


NO 
OFF 

CYLINDRICAL 
.  000( 
50 


YE3 


,  EE4FI1 
.090 'i.l 


Sample  dimensions  : 

A:        .039'+    B:  1B324.0000 


1.9685    D: 


039'+ 
10 


E:  MO 


Spec  i  men  (I 

Max  i  mi  im  1  oad 

Ma::  e  x  t  ens 


I.763 
.047 


Test  end  status  : 

Ma:;  load  point  it 
Ma::  e;:tens  pnt  It 


483 


Number  of  elements 


483 


Specimen  dimensions  : 

A:        .0393701  h :      10324 . 0000O00  C: 


1 . 968504 U    D: 


.039370)     f".  :N0 


Auxiliary    Input    Array    tt    1 

*+<  ****  *     *******    ***■)!*+>-*    *******    **-x-x*>*    *******    *  **  * -x  +  ■«     ^)hh(«.ji*    #*  infjum 

Auxiliary  Input  Array  it  2 

#-X'**  ***   **■**#>'*   *  ***#*■)<   ***#*•*¥   *******   **•*■»»■»+<   >*>****   *-»*.*#.»  *   <HD>*<  »  * 
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GRAPHITE  BUNDLE  TEST  1.5  CM) 


COHOSHES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939<t0 


Test  type:     Yarn/Fiber 
Operator  name:  Mark  Jones 

Sample  Identification:  090902 

Interface  Type:  ^200  Series 

Machine  Parameters  of  test: 

Sample  Rate  (pts/sec):    .33 
Crosshead  Speed  Ua/iin  ):    ,050 


Instron  Corporation 

Series  II  Automated  Materials  Testing  System  v<t.01c 

Test  Date:  September  9,1988 

Sample  Type:   AS-<t 


Humidity  11):    50 
Temperature  (deg.  F):    73 


'BATCH  ID'         090902 
'OIL  OR  NO  OIL'     NO  OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (tex)        2036. 0 
Gauge  length  (en)  5.0000 

Out  of  1  specimens,  0  excluded. 

Displcment  X  Strain     Load 

at  at       at 

Specimen    Maximua  Maximum  Maximum 

Number       (mm)  (X)  .  (KN) 

1      .2500     5.000      .203!) 

Mean:      .2500      5.000      .203<i 

Standard 

Deviation:  - —  

Mean  - 

2.00  ♦  Sdv:  -- --  

Mean  ♦ 

2.00  i  Sdv:  --- -  

Minimum:    .2500      5.000      .203<i 

Maximum:    .2500      5.000      . 203* 
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GRAPHITE  BUNDLE  TEST  (.5  CM) 


C0MPD5I TES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939<i0 


Test  type:     Yarn/Fiber 
Operator  name:  Mark  Jones 

Sample  Identification:  090903 

Interface  Type:  'iP.OO  Series 

Machine  Parameters  of  test: 

Sample  Rate  (pts/sec):    .33 
Crosshead  Speed  (m»/min  ):    .050 


Instron  Corporation 

Series  IK  Automated  Materials  Testing  System  v<t.01c 

Test  Date:  September  9,1983 

Sarple  Type:   AS-1) 


Humidity  I  J  ) :    50 
Temperature  (deg.  F):    73 


'BATCH  ID'         090903 
'OIL  OR  NO  OIL'     NO  OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (ten)        2036.0 
Gauge  length  (cm)  5.0000 

Out  of  1  specimens,  0  excluded. 

Sample  comments:  Maximum  Load  2?. 36  kg  Relaxtaion  Load  to  21.33  Kg 

Displcment  X  Strain  _  Load 

at  at       at 

Specimen    Maximum  Maximum    Ha-isum 

Number       Is*)  (X)       (KN) 

1      .2500      5.000    '  .219(1 

Mean:      .2500      5.000      .21911 

Standard 

Deviation:  

Mean  - 

2.00  «  Sdv:  

Mean  t 

2.00  «  Sdv:  

Minimum:    .2500      5.000      .2191 

Maximum:    .2500     5.000     .219<t 
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6RfiPHiTE  BUNDLE  TEST  (.5  CM) 


COHOSHES  LABORATORY 
NAVAL  FDSTGRADUATE  SCHOOL 
MGNTEREY  CA  93940 


Test  type:     Yarn/Fiber 
Operator  name:  Mark  Jones 

Sample  Identification:  090904 

Interface  Type:  4200  Series 

Machine  Parameters  of  test: 

Sat-ple  Rate  Ipts/sec):    .33 
Crosshead  Speed  (»«i/»in  ):    ,05C 


Instron  Corporation 

Series  II  Automated  Materials  Testing  Syste"1  v1* .01c 

Test  Date:  September  9,1988 

Sasple  Type:   AS-1! 


Humidity  (  X  ):    50 
Teuperature  (deg.  F):   73 


•BATCH  ID'         090904 
'OIL  OR  NO  OIL'     NO  OIL 

Dmensicns : 

Spec.  1 

Lin.  Density  (tex)        2036.0 
Gauge  length  (•■)  5.0000 

Out  of  1  specmens,  0  excluded. 

Saaple  conirents:  Special n  is  slipping  and  test  stopped  at  .46  »« 

Dispkf.ent  X  Strain  Load 

at  at  at 

Spec  isen    Maxiitui  Maxnun  Maxinun 

Nu-ber       Its)  (X)  UN) 

1      .2700      5.400      .2357 

Mean:      .2700      5.400      .2357 

Standard 

Deviation:  

Mean  - 

2.00  i  Sdv:  

Mean  * 

2.00  I  Sdv:  

Mininun:    .2700      5.400      .2357 

Haxi*u«:    .2700      5.400      .2357 
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GRAPHITE  BUNDLE  TEST  (.5  CM) 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939W 


Test  type: 


Yarn/Fiber 


Operator  name:  Mark  Jones 


lnstron  Corporation 

Series  IX  Automated  Materials  Testing  System  v'i.OIc 

Test  Date:  September  9 , 1988 


Sample  Identification:  090905 
Interface  Type:  <t200  Series 
Machine  Parameters  of  test: 
Sample  Rate  (pts/sec) : 
Crosshead  Speed  (»»/«in  ) : 


.050 


Sanple  Type:   AS-1* 


Humidity  IX):    50 
Temperature  (deg.  Fli    73 


'BATCH  ID'         090905 
'OIL  OR  NO  OIL'     NO  OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (tex)        2036.0 
Gauge  length  (ti)  5.0090 

Out  of  1  specimens,  0  excluded. 
Sample  cements:  Sample  slipped 

Displcment  X  Strain  _  Load 

at       at  at 

Specinen    Maximum  Maxinum  Maximum 

Number       (mm)      (XI  UN) 

1      .2900      5.800      .2992 

Mean:      ,2900      5.800      .2992 

Standard 

Deviation:  

Hean  - 

2.00  i  Sdv:  

Mean  ♦ 

2.00  »  Sdv: -- -  

Minimum:    .2900      5.800      .2992 

Maximum:    .2900      5.800      .2992 
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COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  <?39<tO 


GRAPHITE  BUNDLE  TEST  (.5  CM) 


Test  type:     Yarn/Fiber  Instron  Corporation 

Series  II  Automated  Materials  Testing  Systei  v4.0lc 
Operator  na«e:  Mark  Jones  Test  Date:  Septe«ber  9,1998 

Sarple  Identification:  090906  Saiple  Type:   AS-* 

Interface  Type:  4200  Series 

Machine  Paraieters  of  test: 

Saiple  Rate  (pts/sec):    .33  Humdity  I  J  ):    50 

Crosshead  Speed  (m/iin  ):    .050  Tetperature  (deg.  F):    73 

'BATCH  ID'         (WHjl 
'OIL  OR  NO  OIL'     NO  OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (tei)        2034,0 

Gauge  length  lie)  5.0000 

Out  of  I  specirens,  0  excluded. 
Saiple  coeients:  SAMFLE  SLIFFED 

Displcient  K  Strain  ,  Load 

at       at  at 

Specxen    Maxicit  Maxim  Ma«itui 

Nuiber       lis)       (XI  ON) 

I      .2400      4.800      .2664 

Mean:      .2400      4.800      .2664 

Standard 

Deviation:  

Mean  - 

2.00  i  Sdv:  

Mean  ♦ 

2.00  »  Sdv:  

Minnui:    .2400      4.800      .2664 

Mamui:    .2400     4. BOO     .2664 
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GRAPHITE  BUNDLE  TEST  (5  CM) 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939<t0 


Test  type: 


Yarn/Fiber 


Operator  name:  Hark  Jones 

Sapple  Identification:  090907 
Interface  Type:  4200  Series 
Machine  Parameters  of  test: 

Sample  Rate  Ipts/sec):   3.33 
Crosshead  Speed  ( em/ei n  ) : 


lnstron  Corporation 

Series  II  Automated  Materials  Testing  System  v'i.OIc 

Test  Date:  September  9,19B8 

Sample  Type:   A3-<i 


Humidity  I  X  ):    50 
Temperature  (deg.  F):    73 


'BATCH  ID'         090707 
'OIL  OR  NO  OIL'     NO  OIL 

Dimensions: 

Spec.  1 

Lin,  Density  Itex)        2036.0 
Gauge  length  (en)  50.000 

Out  of  1  specimens,  0  excluded. 

Sample  coE*ents:  SAMPLE  SLIFFED  AT  END  OF  TEST  TEST  STOPPED  BY  OPERATOR 

Displcment  X  Strain  ,  Load 

at  at       at 

Specimen    Maximum  Maximum  Ma-mum 

Nurber       (ml  (X)       UN) 

1      .6500      1.300      .2024 

Mean:      .6500      1.300      .2026 

Standard 

Deviation:  

Mean  - 

2.00  »  Sdv:  — -  

Mean  ♦ 

2.00  i  Sdv:  — - - 

Minimum    .6500      1.300      .2026 

Maxieus:    .6500      1.300      .2026 
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GRAPHITE  BUNDLE  TESI  (5  CM) 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93910 


Test  type:     Yarn/Fiber 
Operator  naae:  hark  Jones 

Sai-ple  Identification:  090908 

Interface  Type:  <i200  Series 

Machine  Parateters  of  test: 

Sample  Rate  Ipts/sec):   3.33 
Crosshead  Speed  (»m/ain  ):    .500 


Instron  Corporation 

Series  II  Autonated  Materials  Testing  Systen  v'i.OIc 

Test  Datei  Septenber  9,1988 

5a»p  1  e  Type:   A5-<i 


Humid) t y  (  X  ):    50 
Tetperature  Ideg.  Fl:    73 


'BATCH  ID'         090908 
'OIL  OR  NO  OIL'     NO  OIL 

Dicensions: 

Spec.  I 

Lin.  Density  Hex)        5036. P 
Sauge  length  (in)  50.000 

Out  of  1  specimens,  0  excluded. 

Displcnent  I  Strain  Load 

at  at  ,  at 

Specioen    Maxieut  Maxieuo  Maxinun 

Nueber       («n)  IX)  (KN) 

1      .mo  I.2B0      .2299 

Mean:      .MOO      1.280      .2299 

Standard 

Deviation:  

Mean  - 

2.00  *  Sdv:  

Mean  ♦ 

2.00  *  Sdv:  — - 

Minium:    .MOO      1.280      .2299 

Maxiuuj:    .6W0      1.280      .2299 
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VERSION   4.01c 

Sample  id  :  09090B 
Version  :  k  .01c 
Report  file«  :  12 


Tf?Gt  date? 

Vers  i  on  date 

Operator 


September   9.  19E1D 
06/ 09 /SB 
Mart:  Jones 


X  conversion  :    .03937008 
Y  conversion  :   E.E046S300 


X  A/D  offset 
Y  A/D  off  sot 


,  ("lOOO 


Samp le  rate 

A/D  range 

Cal ib  typo 

Calib  load 

Tempera bure 

Test  type 

B  a  r  t  y  p  p 

Break  ched 

Load  limit 


3.33 
0 
AUTOMATIC 

1 10.231 1 
73 
DEWIER 
AS-4 

1  (!) .  00000 
1  IS. 40450 


E>:  tensome  ber 

Ali  tost  art 

Geometry 

Ca  lib  e>:  tons 

Hum i d i  t  y 

tt  spec ] mens 

Entry  d  i  mens 

Thresh  delay 

E::  tons  limit 


NU 
OFF 

CYLINDRICAL. 
.  0000 


1 

YES 


224B1 
,09843 


B  a m [  1 1 o  c 1 3  mensj o 1 1 s  : 

A:        .0394    B:  1832''+.  0000    C: 


1  .968: 


D: 


.039'*    E:  Ni) 


Spec i  men  tt 

Ma::  i  mum  1  oad 

Ma;:  o::  tot  is 


51  -A'>n 
.044 


Tost  end  status  ! 
Max  load  point  It 
Ma>:  extens  pnt  It 


Number  of  elements  :    450 


Specimen  dimensions  : 
A  :        .03 9 3 7 0  1  B  :  183 S '♦  .  0 0 0 00 0 0  C  : 


1  .  96850'fO  D: 


.0393701  E:ND 


Auxiliary  Input  Array  tt  1 

*•*  *■*  #*■-«   J*******   *.!(.* -tut-.*  *   *  *  *  .*.  .*.  *  *   »)(»innn(   ft.*.*.*.*.*.*   (Humnnt   **».*.»  i,  *   >  *  *  . 

Auxiliary  Input  Array  tt  2 
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GRAPHITE  BUNDLE  TEST  (5  CM) 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939<t0 


Test  type:     Yarn/Fiber 
Operator  nate:  Mark  Jones 

Sample  Identification:  090909 

Interface  Type:  WOO  Series 

Machine  Parareters  of  test: 

Sa»ple  Rate  (pts/sec):       3.33 
Crosshead  Speed  Ue/«in  );         ,500 


Instron  Corporation 

Series  II  Automated  Materials  Testing  System  vMIc 

Test  Date:  September  9,1986 

Sample  Type:   AS-1! 


Humidity  IX):    SO 
Tefperature  (deg.  F):    73 


'BATCH  ID'         090909 
'OIL  OR  NO  OIL'     OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (tex)        2036.0 
Gauge  length  (en)  50.000 

Out  of  1  specimens,  0  excluded. 
Sample  cos'ents:  OIL  SAMPLE 

Displcnent  X  Strain  Load 

at  at  at 

Spec  iffsn    Maxieum  Naxinum  Maxirui 

Number       (mm)  (X)  UN) 

1  .6500  1.300  .237<t 

Mean:      .6500      1.300      .237*1 

Standard 

Deviation:  -- 

Mean  - 

5.00  »  Sdv:  

Mean  f 

E.OO  »  Sdv:  

Minirum:    .6500      1.300      .lllh 

Maximui:    .6500      1.300      .237'. 
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VERSION   '♦.Olc 

Sample  id  :  090909 
Vers  ion  :  *t  .  0  1  c 
Report  filett  :  IS 


Test  date  :  September   9,1980 
Version  date  :  06/09/88 

Operator  :  Marl.  Jones 


X  conversion  : 
Y  conversion  : 


. 03937008 
8.  20'+6S300 


X  A/D  offset 
Y  A/D  offset 


,  (_)IJ<  1(1 

.  (')OOO 


Samp le  rate 

A/D  range 

C  a  1  i  b  type 

Ca lib  load 

Temperature 

Test  type 

Bar  type 

Break  check 

Load  limit 


3.33 
0 
AUTOMATIC 

1 10.831 1 
73 
DENIER 
AS-4 

1 0 . 00000 
1  12.40'-f50 


E::  t  ensomi?  ter 

An  to star t 

Geometry 

Ca  lib  e>:  tens 

Hum  id  i  ty 

tt  spec  i  mens 

Entry  dimens 

Thresh  delay 

E::  tens  limit 


NO 

OFF 

CYLINDRICAL 

.  0000 
50 

1 
YES 

.8S<t01 
.09 8 b3 


Sample    dimensions     : 

A:  .039'+  B:     18324.0000  C; 


1  .9605 


D: 


.  039'+  E:     nn 


Spec  i  men    It 
Ma  :;  i  mi  im     1  oad 

Na"  extens 


53.895 

.  059 


Test  end  status 
Ha::  load  point  II 
Max  extens  pnt  H 


10 


;v,i 
6 1 13 


Number  of  elements  : 


603 


Specimen  dimensions  : 
A:        .0393701  B:  1838'+ .  0000000  C: 


1 .96850^0  D: 


,039370)  E:N0 


Auxiliary  Input  Array  11  1 

#**##■«•*   *<****#*   *  *  *  *  *  *  *-   «.*«*..>.**   ft.*.*.**.**   ^.JH(.JHH(-*   #  *  *  # *  *  *   **  X  *  X  *■*   *  *  *  4,  ■!(  *  * 

Auxiliary  Input  Array  U  8 

****+<* -X  #***■)»•##  *«»»K*  **#*##*  *****##  ###*#*#'  *  *-H  *  *  *  #  Kl'IHH'H   *  »  *  *  *  »  * 
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GRAPHITE  BUNDLE  TEST  0.5  CM  I 


COMPOSITES  LABORATORY 
NAVAL  P05T6RADUATE  SCHOOL 
MONTEREY  CA  93940 


Test  type:     Yarn/Fiber 
Operator  name:  Hark  Jones 

Sasple  Identification:  090910 

Interface  Type:  4200  Series 

Machine  Paraieters  of  test: 

Sample  Rate  (pts/sec):   1,67 
Crosshead  Speed  (am/min  ):    . 250 


Instron  Corporation 

Series  IX  Automated  Materials  Testing  Systet  v<i.01c 

Test  Date:  Septeiber  9,1988 

Sa»ple  Type:   AS-'t 


Huiidity  (  X  ):    50 
Temperature  (deg.  Fl: 


73 


'BATCH  ID'         090910 
'OIL  OR  NO  OIL'     NO  OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (texl        2036.0 
Gauge  length  (ml  55.000 

Out  of  I  specimens,  0  excluded. 

Displcnent  X  Strain  Load 

at        at  ,  at 

Specimen    Maxmu»  Maxims  Maxieu* 

Number       (««)      (X)  UN) 

1      .WOO      1.840      .E174 

Mean:      .WOO      1.840      .S174 

Standard 

Deviation:  

Mean  - 

2.00  I  Sdv:  — - —  

Mean  * 

2.00  »  Sdv:  — ■ 

Minium;    .WOO     I.B40     .E171 

Haxisun:    .WOO      1 .640      .2174 
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GRAPHITE  BUNDLE  TEST  (2.5  CM) 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93W 


Test  type: 


Yarn/Fiber 


Operator  nane:  Hark  Tones 

Sanple  Identification:  100901 
Interface  Type:  'iPOO  Series 
Machine  Paraieters  of  test: 

Saople  Rate  (pts/sec ) :   I 
Crosshead  Speed  ( mm/m 1 n  ): 


.250 


Instron  Corporation 

Series  IX  Autoiated  Materials  Testing  Systei  v<t . 0 1  c 

Test  Date:  Septenber  10,1988 

Sample  Type:   AS-<i 


Humdity  I  X  ):    50 
Temperature  (deg.  Fl:    73 


'BATCH  ID'         100901 
'OIL  OR  NO  OIL'     NO  OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (tex)        2036.0 
Gauge  length  (*m)  £5.000 

Out  of  1  speciiens,  0  excluded. 

Displcnent  X  Strain  Load 

at  at  at 

Speciren    Maxirui  Kaxitua  Maxinui 

N'juber       fin)  (XI  UN) 

1      .«00      1.720      .2130 

Mean:  ,<|300  1.720  .2130 

Standard 

Deviation:     

Mean  - 

2.0"  »  Sdv:  — -     

Mean  * 

2.00  t  Sdv:   

Minnui:         .«00  1.720  .2130 

Maxitun:         .<i300  1.720  .2130 
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VERSION   4.  Olc 

Sample?  id  :  100901 
Vers  1  on  :  4  .  0  1  c 
Repor  t  f  ).  1  e+t  :  1 4 


Test  date 
Version  date 

Operator 


September  10.] 90S 
06/09/88 

Marl:  Jones 


X  conversion 
Y  conversion 


. 03937008 
E.E046E3O0 


X  A/D  offset 
Y  A/D  offset 


.  'X'ui.i 
.  0000 


Gamp  I  e  r  n  k(-: 

A/D  range 

Cal  i  b  type? 

Ca lib  load 

Temperature 

Test  type 

Bar  type 

Break  check 

Load  limit 


1  .67 
0 
AUTOMATIC 

1 10.231 1 
73 
DENIER 
AS-4 

1  0 . 00000 
HE.  404 50 


E::  t  Gnsomn  ter 

Autos tar t 

Geome try 

Cal ib  extens 

Hum  i  ri i ty 

tt  specimens 

Entry  dimens 

Thresh  delay 

E::  tens  limit 


NO 

OFT 

CYLINDRICAL 

.  0000 
50 


YES 


1 . 1E404 
.049 El 


Sample  dimensions 
A :        . 0394    B : 


13324.000 


.9843 


D: 


.0394 


NO 


Spec,  i  men  It 

Ma  ::  \  mum  load 

11a::  extmics 


47 .850 
.  029 


Test  end  status 
Ma::  load  point  It 
Max  e::tens  [int  It 


i  no 
3U5 


Number  of  elements 


303 


Sp ec i men  d  i mens i o ns  : 
A:        .0393701  B:  18324.0000000  C: 


984  2520  D: 


0393701  I":  I  K) 


Aui'i  1  iary  Input  Array  ft  1 

*  *  *  *  *  *•»(   *******   j<  .|H<  *  .*  .)<-  .*   +  *  *.  *  y  «.  *  y  .<  .p  .y.  «.  *  *   *  .»  -x  *.»  **  *  M-¥  *  v.*  *.      +)(fimni   » -K  .  *  *  *  *  Jl 

Auxiliary  Input  Array  tt  2 

*******   *■■*■*•»-  *  **   *******   *******  ********   *  ■¥■**»*  *   **  +  >■*-»*   *»•».•<<->(**   *-»-»->■>-*» 
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GRAFHITE  BUNDLE  TEST  (2.5  CM) 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939<i0 


Test  type: 


Yarn/Fiber 


Operator  name:  Mark  Jones 

Saaple  Identification:  100902 

Interface  Type:  <i200  Series 

Machine  Para«eters  of  test: 

Sample  Rate  (pts/sec):   1.67 
Crosshead  Speed  (nn/cin  ):    .250 


Instron  Corporation 

Series  II  Autoiated  Materials  Testing  System  v<i.0!c 

Test  Date:  Septenber  10,1998 

Sample  Type:   AS-^ 


Huaidity  (  X  ) :    50 
Teeperature  (deg.  F)i    73 


'PATCH  ID'         100902 
'OIL  OR  NO  OIL'     OIL 

Direnstons: 

Spec.  1 

Lin.  Density  (tex)        2036.0 
Gauge  length  (»«)  25.000 

Out  of  1  specimens,  0  excluded. 
Sajple  contents:  SAMPLE  APPEARED  TO  SLIP 

Displcment  X  Strain  Load 

at       at  at 

Specimen    Haxiius  Maxiium  Maxirum 

Nueber       (en)      (X)  •  (KN) 

1      .5000      2.000      .1922 

Mean:      .5000     2.000     .1922 

Standard 

Deviation:  

Mean  - 

2. CO  i   Sdv:  -- -  

Mean  ♦ 

2.00  i  Sdv:  

Minimus:    .5000      2.000      .1922 

Haxinun:    .5000      2.000      .1922 
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GRAPHITE  BUNDLE  TEST  (2.5  CM 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93910 


Test  type: 


Yarn/Fiber 


Operator  nase:  Mark  Jones 

Sample  Identification:  100903 

Interface  Type:  1200  Series 

Machine  Paraoeters  of  test: 

Saeple  Rate  lpts/s?c):   1.67 
Crosshead  Speed  (m/mn  ):    .250 


lnstron  Corporation 

Series  II  Automated  Materials  Testing  System  vl.Olc 

Test  Date:  Septeiber  10,1988 

Saiple  Type:   AS-*t 


Huaidity  I  i   ):    50 
Temperature  Ideg.  F):    73 


'BATCH  ID'         100903 
'OIL  OR  NO  OIL'     OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (tex)        E036.0 
Gauge  length  (id)  25.000 

Out  of  1  specirens,  0  excluded. 
5a»ple  coeients:  SAMFLE  SLIFFED 

Displcirent  %  Strain  Load 

at  at  •     at 

Speciicen          Max  i  nun  Maxuui  Maxuun 

Nusber                 im)  U)  IKN1 

1  .1300  1.720  .2051 

Mean:  .1300  1.720  .2051 

Standard 

Deviation:     

Mean  - 

2.00  f  Sdv:  

Mean  * 

2.00  »  Sdv:  

Nininue:        .<i300  1.720  .2051 

Maxmus:        .1300  1.720  .2051 
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CDHPOSI TES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93940 


GRAPHITE  BUNDLE  TEST  (50  CM) 


Test  type:     Yarn/Fiber  Instron  Corporation 

Series  IX  Autoitated  Materials  Testing  Syste«  v<i.01c 
Operator  na»e:  Mark  Jones  Test  Date:  September  10,1 9B8 

Sample  Identification:  100904  Sarple  Type:   AS-4 

Interface  Type:  4E00  Series 

Machine  Para«eters  of  test: 

Sarple  Rate  (pts/sec):   6.67  Hundity  I  X  ) :    50 

Crosshead  Speed  Imi/iin  ):   5.000  Tenperature  (deg.  F):    73 

'BATCH  ID'         100904 
'OIL  OR  NO  OIL'     NO  OIL 

Dmensions: 

Spec.  1 

Lin.  Density  Itex)        E036.0 
Gauge  length  («»)  500.00 

Out  of  1  specimens,  0  excluded. 

Displcnent  X  Strain  Load 

at        at  at 

Speci«en    Maxitut  Max  i  sum  Max'ifun 

Nuiiber       (sal      IX)  IKN) 

1      3.370      .6740      .1090 

Mean:      3.370      .6740      .1090 

Standard 

Deviation:  

Mean  - 

2.00  ♦  Sdv:  -  - — 

Mean  ♦ 

2,00  ♦  Sdv:  — 

Niniiw:   3.370     .6740     .1090 

Haxinufli:    3.370      .6740      .1090 
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VERSION   4.01c 

Sample  id  :  100904 
Version  :  4.01c 
Report  filett  :  11 


Test  date 
Version  date 

Operator 


September  10,1  90S 

06/09/00 

Mark  Jones 


X  conversion 
Y  conversion 


.03937000 
2.20462300 


X  A/D  offset  : 
Y  A/D  offset  i 


.  0000 


Sample  rate 

A/D  range 

Cal i  b  type 

Calib  load 

Temperature 

Test  type 

Bar  type 

Break  check 

Load  limit 


6.67 
0 
AUTOMATIC 

1 10.531 1 
73 
DEN  I E R 
AS-4 

1 0 . 00000 
1  IS. 4 0450 


E::  tensometer 
Autostar t 

Geometry 
Cal ib  extens 

Hum  id  i  ty 

tt  spec  i  mens 

En t  ry  d i mens 

Thresh  delay 

E::  tens  limit 


NO 

OFF 

CYLINDRICAL 

.  0000 
50 

1 
YES 

.PS'iOl 
.90425 


Sample  dimensions  : 

A:        .0394    B:  10324.0000 


19.6050 


0394 


E :  NO 


Spec  inien  It 

M.i ::  i  mum  1  ond 

11a:;  extens 


24 .472 

.295 


Test  end  status 
N,-i ::  load  point  It 
Ma::  extens  pnt  tt 


10 


271 
601 


Number  of  elements 


601 


Spec i men  d i mens i o ns  : 
A:        .0393701  B:  1 0324 . 0000000  C: 


19.6050400  D: 


.0393701  E:N0 


Auxiliary  Input  Array  it  1 

■>»+-»>*<   4  »  *.  4  *  «  *   -X-J**-***)"   ***  *  *  * *   *******   *****■)»#   ******  -X   *■*■>!*-»*»-   **  *  *  *  4  4 

Auxiliary  Input  Array  It  P. 

*  ¥  -t  I   »  ■»  >   *  *******   4  t>  *  X  *  *  *   *******   *******   *******   *******   ******  *   X  »  »  X  »  x  » 
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GRAPHITE  BUNDLE  TEST  (50  CM) 


COMPOSITES  LABORATORY 
NAVAL  POSTSRADUATE  SCHOOL 
MONTEREY  CA  ?39<i0 


Test  type:     Yarn/Fiber 
Operator  nase:  Mark  Jones 

Saaple  Identification:  100905 

Interface  Type:  <i200  Series 

Machine  Parareters  of  test: 

Sample  Rate  Ipts/sec) :   6.67 
Crosshead  Speed  (*«/«in  ):   5.000 


Instron  Corporation 

Series  IX  Autonated  Materials  Testing  Syste*  v^.Olc 

Test  Date:  Septeuber  10,1988 

Sasple  Type:   AS-1! 


Huiidity  I  X  ):    50 
Temperature  (deg.  P): 


73 


'BATCH  ID'         100905 
'OIL  OR  NO  OIL'     OIL 

Dimensions: 

Spec,  1 

Lin.  Density  (tex)        2036.0 
Gauge  length  Its)  500.00 

Gut  of  1  specimens,  0  excluded. 
Sasple  coc Rents:  OIL  SAMFLE 

Displcment  1  Strain  Load 

at  at  at 

Specimen    Maxirun  Maximum  Maximum 

Number       (in)  (XI  (KNI 

I      3.910      .7820      .1399 

Mean:      3.910      .7820      .1399 

Standard 

Deviation:  

Mean  - 

2.00  »  Sdv: — ~- 

Mean  ♦ 

2.00  «  Sdv:  

Minimum:    3.910      .7820      .1399 

Maxiium:   3.910     .7820     .1399 
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VERSION   4.01c 

Sample  id  :  100905 
Vers i  on  :  4.01c 
Report  filett  :  11 


Tcs  t    date? 

Version    date 

Operator 


September  10, 1980 
06/09/88 

Mart'.  Jones 


X  conversion  : 

Y  conversion  : 

Samp le  rate  : 

A/D  range  : 

Ca 1  lb  type  : 

Calib  load  : 

Temperature  : 

Test  type  : 

Bar  type  : 

Break  check  : 

Load  limit  : 


. 03937008 
E. 80468300 

6.67 
0 
AUTOMATIC 

1 10.231 1 
73 
DENIER 
AS-4 

1 0 . 00000 
1 18.40450 


X  A/D  offset  : 

.  0000 

Y  A/D  offset 

.  0000 

E::  tensometer 

NO 

Au tost sr t 

OFF 

Geometry 

CYLINDRICAL 

Calib  e>:tens 

.  0000 

Hum  i  d i t y 

50 

tt  specimens 

1 

Entry  d l mens 

YES 

1 hresh  del  ay 

.88481 

E::  tens  1  i  in  i  t 

.90485 

Samp  )  e  d  i  mens  i  oris  : 

A:        .039'+    B:  18384.0000 


19.6850    D: 


,0394    E:  NO 


Spec ) men  tt 

Ma::  i  mum  1  oad 

Ma::  extens 


31 .471 
.316 


lest  end  status  : 
Ma::  load  point  1* 
Ma;:  extens  pnt  ti 


10 


314 
644 


Number  o f  elements 


644 


Specimen  dimensions  : 
A  :        . 039370 1  B :  1 0324 . 0000000  C : 


19.6050400  D: 


03937VM  C:NC) 


Auxiliary  Input  Array  tl  1 

»«< I-****   *******   ********   »#*#*#!(   *******   *******   *******   *******   ******* 

Auxiliary  Input  Array  11  2 

»«*>««>'   *******   *******   •»****■«*   *******  **■*  *+*■*'    M  ■»*■¥  **■*■      *  *  *  •*  X  *•  *   ******* 
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GRAPHITE  BUNDLE  TEST  (50  CM) 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939<t0 


Test  type:     Yarn/Fiber 
Operator  name:  Mark  Jones 

Sample  Identification:  100906 

Interface  Type:  WOO  Series 

Machine  Paraeeters  of  test: 

Sanple  Rate  (pts/sec):   6.67 
Crosshead  Speed  U«/«in  ):   5.000 


Instron  Corporation 

Series  II  Automated  Materials  Testing  System  v^.Olc 

Test  Date:  Septeiber  10,1989 

Saiple  Type:   AS-<( 


Humdity  I  X  I:    50 
Tenperature  (deg.  F):    73 


BATCH  ID'         100906 

OIL  OR  NO  OIL'     OIL 

heensions: 

Spec. 

Lin.  Density  (ten) 

5036.0 

Gauge  length  Inn) 

500.00 

Out  of  1  specinens,  0  excluded 


Bisplcjent   X  Strain     Load 
at        at        at 
Specisen    Maxiiuo    Maxi«u«    Maxnui 
Nunber       (in)      (X)       IKN) 

1  3.660      .7320  .1359 

Mean:  3.660      .7320  .1329 

Standard 

Deviation:  -  

Mean  - 

2.00  i  Sdv: — - - - 

Mean  * 

2.00  *  Sdv:  —  - 

Minimi*!  3.660     .7320  .1329 

Maxieun:  3.660      .7320  .1329 
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GRAPHITE  BUNDLE  TEST  (25.0  CHI 


COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939<i0 


Test  type:     Yarn/Fiber 
Operator  na»e:  Mark  Jones 

Saaple  Identification:  10090? 

Interface  Type:  <i200  Series 

Machine  Paraneters  of  test: 

Sarple  Rate  Ipts/sec):   5.00 
Crosshead  Speed  (««/»in  ):   2.500 


Instron  Corporation 

Series  IX  Autoaated  Materials  Testing  System  v<i.01c 

Test  Date:  Septenber  10,1988 

Sample  Type:   AS-** 


Humidity  I  X  I:    50 
Temper  a ture  Ideg.  PI:    73 


'PATCH  ID'         100907 
'OIL  OR  NO  OIL'     OIL 

Di»ens>ons: 

Spec.  1 

Lin.  Density  (tex)        2036.0 
Gauge  length  (ml  250.00 

Out  of  1  specimens,  0  excluded. 

Displceent  X  Strain  Load 

at  at  at 

Specimen    Maxi«u«  Maxirun  Maximum 

Number       (mm)  (X)  (KN) 

1      2.370      .WO  .1716 

Mean:      2.370      .9<t80      .1718 

Standard 

Deviation:  

Mean  - 

2.00  J  Sdv:  - - 

Me^n  ♦ 

2.00  •  Sii:   - 

Minimum:    2.370      .MBO      .1718 

Maximum:    2.370      .9*180      .1718 
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VERSION   ^. Olc 
Sample  id 
Vers  ion 
Repor t  f  i 1 e# 

X  conversion 
Y  conversion 


1 00907 
<+ .  0 1  c 

15 

.03937008 
2.20462300 


Test  date  :  September  Id, 1980 
Version  date  :  06/09/88 

Operator  :  Mart  Jones 


X  A/D  offset 
Y  A/D  offset 


.  OU U(.) 

.  0000 


Samp ] e  rate 

A/D  range 

Ca lib  type 

Ca lib  1 oad 

Temperature 

Test  type 

Bar  type 

Br ea I'  check: 

L  o  a  d  .1  i  m  i  t 


AUTOMATIC 

1 10.2311 
73 
1)  F.N  IK  ft 

1 0 . 00000 

IIS.  '+0't50 


E>:  tensome  t  er 

Autos tar  t 

Geometry 

Ca  lib  e>:  tens 

Hum  i  d  i  ty 

it    spec  i  mens 

Entry  d  i  mens 

Thresh  delay 

E>:  tens  limit 


NO 

OFF 

CYLINDRICAL 

.  0000 
50 

1 
YES 

.22  '+81 
.'+9213 


Sample    rJimr'vinions     i 

A:  .  039't  B:     1032'+.  0000 


9.0'+RI 


Ds 


Spec  i  men    tt     :  1 

M a : :  i  m u m     1  o  nil     :  3 8  .  5 7 3 

Ma::    exterts     :  .179 

Number     of    elements     :  5't? 

Specimen    dimensions     : 


A: 


.0393701  B:  1 032'+ .  0000000  C: 


Test  end  status  : 
Max  load  point  0 
Ma::  ex  tens  put  0 


9.8425200  D: 


039'+ 
10 


E:  NO 


283 
5'+7 


.0393701  E:N0 


Auxiliary  Input  Array  tt  1 

*******  *******  *******    *******  *******  *******  *******  *******  *  *  *  *  *  *  * 

Auxiliary  Input  Array  tt  2 

*******  *******  *******  *******  ********  *******  *  *  *  *  *  *  *  *******  ******* 
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GRAPHITE  BUNDLE  TEST  (25.0  CM) 


COMPOSITES  LABORATORY 
NAVAL  P0ST6RADUATE  SCHOOL 
MONTEREY  CA  939W 


Test  type:     Yarn/Fiber 
Operator  name:  joe 

Sasple  Identification:  100908 

Interface  Type:  <i200  Series 

Machine  Paraseters  of  test: 

Sanple  Rate  (pts/sec):   5.00 
Crosshead  Speed  (nt/nn  ):   E.500 


Instron  Corporation 

Series  IK  Automated  Materials  Testing  System  v4.0lc 

Test  Date:  September  10.19BB 

Saiple  Type:   AS-'i 


Hundity  (  I  ):        50 
Te»perature  (deg.  F):        73 


'BATCH  ID' 

100903 

'OIL  OR  KO  OIL' 

NO  OIL 

Dimensions: 

Spec. 

Lin.  Density 

(tex) 

2036.0 

Gauge  length 

dim) 

250.00 

Out  of  1  specimens,  0  excluded 


Displc»ent  X  Strain  Load 

at  at  at 

SpecUen  Maxuun  Maxi»u»  -Maxirut 

Nusber  (ma)  IX)  (KN) 

1  2.370  .9560  .169d 

Mean:  2.390  .9560  .169<( 

Standard 

Deviation:  

Mean  - 

2.00  ♦  Sdv:  

Mean  * 

2.00  I  Sdv:  

MiniBu*:  2.390  .9560  ,169't 

Maxnui:  2.390  .9560  .I6?<i 
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VERSION   4.01c 

Sample  id  :  100900 
Version  :  4.01c 
Report  f  i 1 ett  :  15 


Test  date 
Version  date 

Operator 


September  10. 1980 

06/09/00 

joe 


X  conversion 
Y  conversion 

Sample  rate 

A/D  range 

Ca lib  type 

Calib  load 

Temperature 

Test  type 

Bar  type 

Break  check 

Load  limit 


,03937000 
,204 62300 


AUTOMATIC 

1 10.231 1 
73 
DENIER 
AS-4 

1 0 . 00000 
1 12.40450 


X  A/D  offset  : 

.  0000 

Y  A/D  offset  : 

.  0000 

E>:  tensometer 

NO 

Autos  tar  t 

OFF 

Geometry 

CYLINDRICAL 

Cf\  lib  e::  tens 

.  0000 

Hum  i  d  i  ty 

50 

ft  spec l mens 

1 

Entry  d  i  mens 

YES 

Thresh  delay 

.22481 

E::  tens  limit 

.49213 

Sample  dimensions  : 

A:       .  0394    B:  18324  .0000 


9. 8425    D: 


.0394 
1  0 


F. :  NO 


Spec  i  men  -It 

Mas:  i mum  )  oad 

Max  ex tens 


30.01 1 
.  174 


Test  e nd  statu  s 
Max  load  point  0 
ha:;  ex  tens  pnt  tt 


207 

530 


Number  of  elements 


530 


S p  e c  i  me n  cl  i  men s  i  o  n s  : 
A :        .0 3 9 3701  D :  10 3 24 . 0 0 0 0 0 0 0  C : 


9.0425200  D: 


,0393701  E:N0 


Auxiliary  Input  Array  tt  1 

*#»#**«  *******  *******  *******  *******  *******  *******  *******  ******* 

Auxiliary  Input  Array  It  2 

*  x  x  -x  -x  *  x  *  *  -x  *:  -x  »  *  ******>  -x  *  x  -x  x  *  -x  *  *  x  *  *  *  x  *  x  *  *  *  4  *  *  *  *  *  x  x  •*  mm*       .10 
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